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Abstract: We develop techniques to compute multi-instanton corrections to the 1/N expansion 
in matrix models described by orthogonal polynomials. These techniques are based on finding 
trans-series solutions, i.e. formal solutions with exponentially small corrections, to the recursion 
relations characterizing the free energy. We illustrate this method in the Hermitian, quartic 
matrix model, and we provide a detailed description of the instanton corrections in the Gross- 
Witten-Wadia (GWW) unitary matrix model. Moreover, we use Borel resummation techniques 
and results from the theory of resurgent functions to relate the formal multi-instanton series to 
the nonperturbative definition of the matrix model. We study this relation in the case of the 
GWW model and its double-scaling limit, providing in this way a nice illustration of various 
mechanisms connecting the resummation of perturbative series to nonperturbative results, like 
the cancellation of nonperturbative ambiguities. Finally, we argue that trans-series solutions 
are also relevant in the context of topological string theory. In particular, we point out that in 
topological string models with both a matrix model and a large N gauge theory description, the 
nonperturbative, holographic definition involves a sum over the multi-instanton sectors of the 
matrix model. 
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1. Introduction 

The study of nonperturbative effects in the 1/N expansion of matrix models is of great theo- 
retical and practical importance. First of all, they provide a toy model for the nonperturbative 
aspects of the 1/N expansion in more complicated theories. Since matrix models are able to 
describe both noncritical string theories as well as topological string theories on certain back- 
grounds, they provide a reliable arena for nonperturbative computations in string theory. Finally, 
the critical points of matrix models offer a surprising catalogue of scaling behaviors which are 



relevant in many physical and mathematical problems (see [26] for a recent overview), and their 
nonperturbative aspects should be also important in that context. 

Nonperturbative effects in matrix models were identified long ago in terms of eigenvalue 
tunneling |77], E3|, and they were subsequently analyzed in the so-called double-scaling limit 
[46, 30, 12 1 that describes conformal field theories coupled to gravity (see, for example, [36, 23, 
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24, |52j). Surprisingly, much less is known about these effects away from the critical point. In [66] 
various general results were obtained for one-instanton amplitudes in one-cut matrix models, and 
it was explicitly shown in many examples that these instantons govern the large order behavior 
of the 1/N expansion. A very important motivation for this study was the connection between 



matrix models and topological string theories discovered in |2S|] and extended in [65, 10 1 to toric 
backgrounds. In some simple cases, the methods of J66[ ] could be applied in order to obtain 
nonperturbative corrections to the total free energy of topological string theory. 

Another way to compute instanton effects in a matrix model is to consider general, unstable 
multi-cut configurations. This was pointed out in ||, and developed more recently in ]33| in the 
Hermitian case to obtain a formal, universal expansion for the matrix model partition function 
which can be regarded as a multi-instanton expansion. Both |3(| and |33|] are based on the 
geometric description of the 1/N expansion in terms of a spectral curve. This formalism is very 
powerful since it gives universal expressions, but it is not so easy to implement in practice. In 
some cases it is much more convenient to compute the relevant amplitudes in the 1/N expansion 
by using the method of orthogonal polynomials, which was introduced and developed in fj], [g and 
played a key role in the analysis of the double-scaling limit. For example, if one has to compute 
the free energies F g of the quartic matrix model at high genera, the method of orthogonal 
polynomials will be much more efficient than methods based on the spectral curve. One could 
suspect that the same considerations apply to the calculation of nonperturbative effects. 

One of the purposes of this paper is in fact to develop techniques to compute multi-instanton 
effects in matrix models by using the formalism of orthogonal polynomials. The principle behind 
these techniques is rather simple, and it can be easily motivated by considering the double- 
scaling limit of the matrix model. In this limit, the recursion relations of orthogonal polynomials 
lead to a differential equation (usually called the string equation) for the specific heat. An 
asymptotic series solution of this equation gives then the perturbative free energy. However, 
one can also compute multi-instanton amplitudes from the string equation by considering a 
trans-series solution to the differential equation, i.e. a solution involving exponentially small 
corrections to the asymptotics. The resulting amplitudes are double-scaling limits of the full 
multi-instanton amplitudes off-criticality. In order to obtain these, one notices that the string 
equation is obtained from a difference equation, and as we will show in this paper, the trans- 
series solution to this difference equation gives a systematic way to compute the nonperturbative, 
multi-instanton effects of the full matrix model. 

As a first illustration of this method, we revisit a canonical example, namely the quartic 
matrix model, and we recover and extend the results of [66] for the one-instanton amplitude. 
One of the advantages of the techniques we develop here is that it they can be also applied to 
unitary matrix models, where the geometric techniques based on spectral curves have not been 
developed, and we present a detailed analysis of the instanton corrections in the simplest unitary 
model, namely the Gross-Witten-Wadia model [47, p0|] . This model is very interesting in many 
respects: it is the model underlying two-dimensional Yang-Mills theory on the lattice, it has a 
double-scaling limit |0] which describes the simplest minimal superstring theory [59, 76], and 
it plays a role in the description of Yang-Mills theory on S 3 x S 1 (see for a recent study and 
references to previous work). 

The techniques we develop give the multi-instanton corrections as formal power series in 
two small parameters, namely 1/N and e~ N (this is indeed typical of a quantum-mechanical 
computation involving instanton effects, like the WKB method). Once one has determined these 
formal, nonperturbative corrections, a natural question is: how is the resulting series related 
to the nonperturbative definition of the model? The second purpose of this paper is to clarify 
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this relation. It turns out that, in the cases we will be interested in, the trans-series can be 
resummed with Borel resummation techniques, and one obtains in the end a one-parameter 
family of functions which can be regarded as nonperturbative completions of the theory. The 
fact that formal solutions to the relevant equations come in families is very well known in the case 
of double-scaled matrix models, but we use results in the theory of exponential asymptotics and 
the theory of resurgent functions to construct actual (convergent) solutions both at criticality 
and off-criticality. 

We illustrate these ideas in the case of the GWW model and its double-scaling limit. The 1/N 
expansion of this model, in the weak coupling regime, has the properties typical of perturbative 
expansions in realistic quantum field theories: the expansion is not Borel summable, yet the model 
has a unique nonperturbative definition in terms of the original unitary integral. Therefore, this 
is an excellent laboratory to explore the analyticity properties of the 1/N expansion. We show 
that the 1/N series can be resummed in such a way that the ambiguities coming from the Borel 
resummation cancel against nonperturbative instanton effects, providing a nice illustration of the 
cancellation mechanism much discussed in renormalon physics [22, E3] and in some quantum- 



mechanical problems |83|]. In our case, as in U§3 1, this cancellation is a consequence of the 
"resurgence" properties of trans-series expansions [[H], 16]. The resummation process gives a 
one-parameter family of solutions which include in a crucial way multi- instanton corrections. 
We show that, for a particular (and rather natural) choice of the parameter, this solution is the 
semiclassical expansion of the true nonperturbative answer, and multi-instanton corrections are 
crucial in order to reproduce the exact, nonperturbative value of the different physical quantities. 

Our results for the double-scaling limit of the GWW model have a clear interpretation 
in terms of minimal superstring theory. As in the case of bosonic minimal string theories, 
the minimal superstring has ZZ branes which should correspond to eigenvalue instantons. We 
interpret the sum over multi-instantons that arises naturally in the trans-series solution as a sum 
over ZZ brane backgrounds, i.e. over sub-leading saddles of the theory. It was pointed out in [61] 
that the exact, nonperturbative answer for the free energy in noncritical string theories should 
include indeed a sum over all of these backgrounds. This is precisely what we obtain here. The 
general lesson of our analysis is that, to make sense of this sum, one has to be careful about 
various subtleties. These include Borel resummation of the various asymptotic expansions, the 
choice of nonperturbative parameter, the cancellation of non-perturbative ambiguities, and the 
reality conditions of the final solution. 

As we mentioned before, an important motivation for the study of nonperturbative effects 
in matrix models is the connection to topological string theory. It is natural to conjecture that 
the full topological string theory partition function is a trans-series expansion which includes 
nonperturbative multi-instanton effects. This was already suggested in |3l| following the analogy 
with noncritical strings, and made more concrete in [pq ]. In some simple topological string 
theories with a dual matrix model description, the one- instanton sector was studied in |6(| by 
using their geometric description in terms of spectral curves/mirror symmetry. In this paper we 
give two more pieces of evidence for this conjecture. We first revisit one of the examples of [66], 
namely Hurwitz theory, which can be regarded as a toy model of topological string theory. The 
free energy of this theory is described by a difference equation of the Toda type [70|, therefore 
one can find a trans-series solution of this equation and we verify that this solution reproduces 
the one-instanton effects computed in [66]. Second, we point out that in more complicated 
topological string models, with both a matrix model and a holographic, large N gauge theory 
description, the nonperturbative, holographic definition tells us that the full partition function 
involves indeed a sum over instanton sectors in the matrix model. It has been pointed out by 
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Eynard [^] that such sums should be background independent, since all the backgrounds are 
summed over, and this opens the possibility that a proper understanding of such trans-series 
solutions will lead to background independent topological string models. 

The organization of this paper is as follows. In section 2 we review some basic results on 
matrix models in the formalism of orthogonal polynomials. In section 3 we explain how to 
obtain multi-instanton amplitudes in matrix models by constructing trans-series solutions to the 
relevant difference equations. We illustrate the method for the quartic matrix model. Section 4 
is devoted to a detailed analysis of instanton effects in the GWW unitary matrix model, focusing 
on the weakly coupled phase. In section 5 we address the issue of how to relate the formal multi- 
instanton expansions to the nonperturbative definition of the matrix model (when available). 
We introduce ideas and techniques from Borel resummation and the theory of resurgence, and 
we perform a very detailed analysis of this issue in the GWW model and its double-scaling 
limit, described by Painleve II. Section 6 uses the framework developed in this paper to analyze 



topological string theory, building on the results of [56|. Finally, section 7 states our conclusions 
as well as some open problems. 

2. Matrix models and orthogonal polynomials 

In this section we will review some elementary aspects of the method of orthogonal polynomials 
as applied to the calculation of the 1 /N expansion of a matrix model. Many of the key formulae 
that we will need are common to both Hermitian and unitary matrix models, but for concreteness 
we will start discussing the Hermitian case. The method of orthogonal polynomials for Hermitian 
matrix models was discovered in [[?], ||], and useful reviews can be found in (2^, 64]. We will follow 
the conventions of this last reference. 

We will consider gauged, Hermitian matrix models defined by the partition function 

where V{M) is the potential. A standard argument reduces this integral to an integral over 
eigenvalues 

l/n^A 2 (A)e-^».^> (2.2) 

' J i = l ' 

where 



A(A) = \{{X l - A,) (2.3) 

i<j 

is the Vandermonde determinant. If we regard 

d/i = e -^(A)dA 

as a measure in R, one can introduce orthogonal polynomials p n (A) defined by 

dnp n (\)p m (\) = h n 5 nm , n > 0, (2.5) 



where p n {\) are normalized by requiring the behavior p„(A) = A n + • • • . One then easily finds, 

N-l N 

Z=X\h l = h»\{r?-\ (2.6) 



i=0 i=l 
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where we have introduced the coefficients 

r k = ~r~~) k > 1 (2.7) 

which appear in the recursion relations for the p n (A), 

(A + S n )p n (X) = p n+ i(X) + 

'nPn—1 

(A). (2.8) 
It will be useful to normalize the results by considering the Gaussian matrix model, 

1 f dMe~^ T ™\ (2.9) 



J G 



vol(U(N)) 

i.e. we will be interested in computing the normalized free energy 

F = log Z — log Zq. (2.10) 
This free energy has an asymptotic expansion around g s = of the form 

oo 

F(t,g s ) = ^2F g (t)g 2 r 2 (2.11) 

9=0 

where t is the 't Hooft parameter 

t = g s N. (2.12) 



Since we keep t fixed, ( 2.11 ) is also a large N expansion in powers of 1/N 2 . The standard 
procedure to compute this asymptotic expansion by using orthogonal polynomials goes as follows. 
We have an exact formula for finite N, 

9 t 2 , h Q t 2 v^v { k\ r k , 



N h$ Nfri\ Nj "kg, 

where is the coefficient ho for the Gaussian model. In order to proceed, we introduce a 
continuous variable as N — > oo, 

g s k -> z, < z < t, (2.14) 
and we assume that in this continuum, N — > oo limit, r k becomes a function of z and g s , 

r k ^R(z,g s ) (2.15) 

It will be useful to consider the function 

R{z,g s ) 

a(z,g a ) = (2.16) 

z 

which can be regarded as the continuum limit of r k /{kg s ). It is easy to see that, for polynomial 
potentials of the form 

V(M) = \m 2 + ■■■ (2.17) 
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one has ~ kg s + • • • , therefore the function log(rf./(kg s )) is regular at k = and we can use 
the standard Euler-Maclaurin summation formula to evaluate ( 2.13| ), One then obtains |<], 0]: 



>F= / dz(t-z)\ogZ(z) + V<?. 
/i 



2p B 2p d 2 ^ 1 



(2p)! dz 2 ?- 1 



(t-z) logZ(z,g s ) 



z=t 



z=0 



(2.18) 



+ 



tg s 



^0 



We will rephrase ( p,18| ) in a more convenient way. A small calculation shows that 



,d 2 F 
dt 2 



We now use the fact that 



B 



2p 



d 2 P 



z 2 csch 2 (z) 



o2fc D 



k=l 



{2k){2k - 2)!' 



to write the above equation as 



4 sinh 2 



2 dt) 



F(t) =log! 



The first member can be written as a difference operator, therefore 

F(t + ffs ) + F(t - g s ) - 2F(t) = log H, 

or equivalently, 



exp 



F(t + g s ) + F(t - g s ) - 2F(t)\ = S. 
A shorter way to derive this equation is simply to start from the identity 

Zn+iZn-i 



Z N 



rN, 



(2.19) 
(2.20) 

(2.21) 

(2.22) 
(2.23) 

(2.24) 



where is the partition function ( [2.6| ) at rank N, and consider its continuum limit. Notice 
that, written in the form (2.23), the equation determining the free energy involves the standard 
difference operator of the Toda lattice. This is related to the fact that the free energy of a 
polynomial matrix model is a solution to the Toda hierarchy fjM. 

In order to compute the g s expansion of the free energy (|2.1l[) , one finds first an expansion 
for R(z,g s ) of the form 



R^(z,g s ) = Y j gl s R Qt2s (z)- 



(2.25) 



s=0 



Once this expansion is plugged in E(z,g s ) and then in ( |2.18| ), the expansion ([ETll) follows. In 
order to obtain ( [2.25 ) one has to use the so-called pre-string equation. This is a difference equation 
for R(z,g s ) which can be derived as the continuum limit of the recursion relations obeyed by the 
coefficients (2.7). The pre-string equation can be explicitly written for any polynomial potential 

ic 

1 



28]. For example, in the case of the quartic matrix model with potential 

A 



V(M) 



-M 



48 



-AF 1 



(2.26) 
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the difference equation for R(z,g s ) reads as 

R(z, g s ){l- ^(R(z, g s ) + R(z + g s ,g s ) + R(z - 9s, 9s)} = z. (2.27) 

This type of difference equations have a solution of the form ( |2,25| ), and they determine Rq s (z) in 
terms of the Ro iS >(z), s' < s. When this solution is plugged in fl2.22|) , one obtains the perturbative 
expansion of the total free energy in powers of g s , which is the standard 1/N expansion of the 
matrix model |7], |8| . 

The formalism of orthogonal polynomials for unitary matrix models is very similar (see, for 
example, |ffl| ). We consider unitary matrix models of the form 

Z= ! dUe^ V{U) , (2.28) 



where U is a unitary matrix and the potential V{U) has the structure 

V(U) = J2(9ltrU l + g l txU^, 9l = - i 7l ). (2.29) 

We can write the partition function in terms of the eigenvalues of U, (pi 6 [— 7r, 7t]: 

Z = I JJd^n4sin 2 (^^)e^^ V ^, (2.30) 

J i i<j ^ ' 

v((i))=j2^j cosi( i ) +j sini ^y ( 2 - 31 ) 



where 



If we introduce the measure 

the orthogonal, monic polynomials 
satisfy 



1 dz ivy . . 

dfi = — — e^ (2.32) 
Ztti z 

Pn {z) =/ + •■• (2.33) 

dfj, p n (z)p m (z' 1 ) = h n 5 nm . (2.34) 

as well as the the recursion relation 

p n+ i(z) = zp n (z) + f n z n p n (z' 1 ), (2.35) 

and one easily shows that 

= l-fl (2.36) 



K 

As in the Hermitian case, we introduce the quantities 

r n = (2-37) 
fln-l 

In terms of these, the partition function of the unitary matrix model is given again by the formula 
(|2.6| ). Normalizing by the Hermitian, Gaussian matrix model, and introducing the continuum 
limit Q2-14 ) for the r n , we find the same formalism describing both the unitary and the Hermitian 
matrix model. In both cases the key ingredient is to derive explicit expressions for the function 
R(z,g s ). In the unitary case these are also obtained by solving a difference equation, which 
will depend on the particular model one is considering. A particularly important example, the 
Gross- Witten-Wadia (GWW) model |47|, |o|, will be studied in section 4. 
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3. Multi-instanton corrections 



In this section we develop techniques to compute multi-instanton corrections in matrix models 
with the method of orthogonal polynomials. We first describe the general structure of the method, 
and then, as an example, we analyze in some detail the Hermitian matrix model with an even, 
quartic potential. 

3.1 Multi-instantons and matrix models 

Before presenting our method we will provide a general framework for instanton calculus in 
matrix models, referring to for example [^, ^4|, |5^, ||, ^] for more details. In this section we 
will consider Hermitian matrix models with polynomial potentials. 

The expansion of the matrix model partition function Z in even powers of g s is in fact an 
asymptotic expansion of Z around a saddle point of the matrix integral. These saddle points 
are characterized by a distribution of matrix eigenvalues p(A). For example, in the so-called 
one-cut case, all the eigenvalues sit on the same interval, which is located around a minimum 
of the potential. Let us assume that the potential of the matrix model has d different extrema 
xi, ■ ■ ■ ,Xd- Then, the most general saddle-point is a configuration in which the N eigenvalues 
split into d sets of Nk eigenvalues, k = 1, • • ■ ,d. Let us denote each of these d sets by 

R ( fk=l,..,7V fe , k = l,---,d. (3.1) 

The eigenvalues in the A;-th set sit in an interval or arc Ik around the k-th extremum. Along this 
interval, the effective potential 

V eS (X) = V(X) - t J dA'p(A') log | A - A'| (3.2) 

is constant. It is possible to choose d integration contours Ck in the complex plane, k = 1, • • • ,d, 
going to infinity in directions where the integrand decays exponentially, and in such a way that 
each of them passes through exactly one of the d critical points (see for example The 
resulting matrix integral is convergent and can be written as 

Z{N 1 ,..,N d ) = —^—( -../ l[pA\X)e-^ V ^. (3.3) 
iVi! • • • N d \ J X W 6Cl J x w eCd JLX 2vr 

Of course, when the integrand is written out in detail, it splits into d sets of eigenvalues which 
interact among them through the Vandermonde determinant (see for example [|6C[| ). If one 
now regards (|3.3| ) as the matrix integral in a topological sector characterized by the numbers 



Ni, ■ ■ ■ , Nd, it is natural to consider the general partition function [g^, 24, ^, 33 1 

Z= Ci 1 --<i d Z(Nir-- M. (3.4) 

Ni+-+N d =N 

The coefficients Ck can be regarded as 9 parameters which lead to different 9 vacua ID 1 . Notice 
that one can fix the overall normalization by setting one of the Cfc's to 1, for example. The sum 
( |3.4| ) can be also regarded as a matrix integral where the iV eigenvalues are integrated along the 
contour 

d 

C = Y,CkCk, (3.5) 

k=l 



1 This type of structure has been argued to be relevant as well for general QFT path integrals [|42|, |49| 
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Figure 1: An ^-instanton configuration in a matrix model with a cubic potential. N — I eigenvalues sit 
at the minimum, while I eigenvalues sit at the maximum. 



therefore the 9 parameters give the relative weight of the different contours [£4|, 33]. 

Among all configurations characterized by the fillings (Ni,--- , Nd), the most stable one 
occurs when all the eigenvalues sit at the minimum of the potential, which we will take to 
be x\ (we assume for simplicity that this minimum is unique). We will regard the resulting 
configuration 

( AT, 0, ••-,()) (3.6) 

as the reference configuration for the system. This corresponds to a one-cut solution of the 
matrix model. It is then easy to see that the other terms in the sum ( |3,4D , with general filling 
numbers Ni, are exponentially suppressed with respect to the reference configuration ( ]3.6| ), with 
a weight of the form 

ex P (--y)JV < (V r e ff (x i )-Kff(xi))) (3.7) 

and can then be regarded as instanton configurations. 

An example of such a situation is the cubic matrix model, where the potential has two 
different extrema (a minimum and a maximum). In the reference configuration all the eigenvalues 
of the matrix sit near the minimum. This is the standard one-cut solution described by the 
method of orthogonal polynomials, and gives Z(N,0). The configuration in which I eigenvalues 
sit near the maximum is an ^-instanton of the cubic matrix model, and it gives the partition 
function Z(N — £,£) (see Fig. [l]). After summing over all topological sectors, and fixing the 
normalization by setting = 1, £2 = C> we obtain the partition function 

N 

Z = Y J C e Z(N-e,£), (3.8) 
1=0 

which in the large N limit becomes 

00 

Z = Y J <*Z®{t,g a ). (3.9) 

£=0 

The method of orthogonal polynomials, as presented in 0, ||] and summarized in the previous 
section, gives tools to compute the asymptotic g s expansion of the one-cut answer Z^°\t,g s ), and 
it is natural to ask if we can use it to calculate the instanton corrections to Z^°'{t, g s ), i.e. the 
partition functions Z^'(t,g s ), I 7^ 0. Each of these should have an asymptotic expansion in g s , 
which corresponds physically to the perturbative g s expansion around the instanton configuration. 
In this paper we will develop techniques to do that. We will focus on the situation that we 
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have just considered, namely, a matrix integral where the reference configuration is the one-cut 
solution, and the instantons are obtained by moving a small number of eigenvalues from the 
minimum to another saddle point. 

It is worth mentioning that (|3.3| ) is nothing but a partition function for a multicut matrix 
model, and one should be able to evaluate a generic instanton configuration as a particular 
case of the multicut theory. This was pointed out in || and developed recently in more detail 
in 33], where general expressions were obtained for the formal expansion of (3.3) and ( p.4| ) in 



the Hermitian case, and in the framework of the full saddle-point 1/N expansion presented in 
p5[| . Here we are interested in computing this formal expansion in concrete models described by 
orthogonal polynomials, and to work out the precise relation between the formal expansion and 
the original integral. Some aspects of the relation between multi-instantons and multicut models 
will be studied in J67|]. 

3.2 Multi-instantons and trans-series 

In order to understand our approach to the calculation of nonperturbative effects in matrix 
models off-criticality, it is useful to look first at this problem in the double scaling limit (see 
p8[| for a review of the double-scaling limit of matrix models). For simplicity we will consider 
the example of pure 2d gravity, which can be obtained for example from the quartic or the 
cubic Hermitian matrix models. In this limit, the total free energy F^ s (k), as a function of the 
cosmological constant k, is described by a set of two equations. The first one relates to the 
specific heat u, and reads 

= -u(k). (3.10) 

This is in fact the double-scaling limit of fl2.22| ) . The second equation is a differential equation for 
u. For 2d gravity one obtains Painleve I (we use the normalization appropriated for the quartic 
matrix model) 

--u" + u 2 = k. (3.11) 
3 

The asymptotic solution of this equation which goes like u ~ \fk as k — > oo, 

oo 

(«) = Vk^uo^k 59/2 (3.12) 

3=0 

describes the perturbative free energy. However, one can find a one-parameter family of solutions 
to ( |3.11 ) which includes exponentially suppressed terms as k — > 00: 



u(«) =^CV £ )(k) = ^£Ck-¥ e" MK5/4 eW( K ), (3.13) 
1=0 e=o 

where C is a parameter, the constant A has the value 

A^-f (3.14) 

and 



00 

£ W( K ) = J>,, n+lK - 5n / 4 (3.15) 

n=0 
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are asymptotic series. Since we have introduced an arbitrary constant C in ( |3.13 ), we can 
normalize the solution such that u\ t o = 1. 

These types of solutions to differential equations are called trans-series, and are the central 
object in the theory of exponential asymptotics. By plugging ( |3.13| ) in ( |3.10| ), we obtain a similar 
trans-series expansion for the free energy, 

oo oo 

F ds („) = £ C^j?(K) = £ C^-Se^VSC*), (3-16) 
1=0 e=o 

where Fjj? is the perturbative free energy. The exponentially suppressed corrections in the trans- 
series expansion can not be seen in a perturbative expansion around k = oo and their origin is a 
nonperturbative effect. 



Several remarks can be made concerning these solutions. The first one is that ( 3.13 ) is an 
expansion in two small parameters, namely 

X = K -V\ £ = K -§ e -^ 5/4 . (3.17) 

Usually one first expands in £ in order to extract the i-th. term u™' in the trans-series, and then 
one expands in x in order to obtain the asymptotic expansion of this term, but in some cases it 
is useful to first expand in x [fi9| . The second remark is that ( |3,13j ) is a one-parameter family, 
parametrized by C. This corresponds to the nonperturbative ambiguity plaguing these problems. 
The third remark is that, as for any asymptotic expansion, the trans-series solution is only valid 
in a sector of the complex plane, and as we go from one sector to another and cross a Stokes 
line the asymptotics will change. The Stokes line in the example above is the positive real axis 
arg(ft) = 0. However, the difference between the two asymptotic solutions as we cross a Stokes 



line will be a shift in the parameter C appearing in (3.13) 



C -> C + S, (3.18) 

where S is sometimes called the Stokes multiplier. There are many ways to obtain S in the case 
of Painleve I. One can for example deduce it from a matrix model calculation, as first done by 
David in \24\, or one can derive it rigorously in the framework of isomonodromy deformations, 



see [40 1 and the comprehensive book 

What is the nonperturbative origin of the exponentially suppressed terms in the trans-series 
expansion? In the case of the expansion describing the double-scaling limit of matrix models, 
these effects are due to multi-instantons. The ^-instanton correction F^(k), which is obtained 
from the full trans-series solution for u(n), can be computed by taking the double-scaling limit of 
an £-instanton configuration of the appropriate matrix model partition function. For Painleve I, 
one can take for example the cubic matrix model, and the free energy F^(k) can be computed 
as the double scaling limit of log Z^\t, g s ) appearing in ( |3.9[ ). For I = 1, the instanton origin of 
the exponentially suppressed corrections to Painleve I has been verified by a direct calculation 



in [23, 24, 52, pq |, in the context of the saddle-point solution to Hermitian matrix models. 

Of course, it is much more efficient to compute multi-instanton effects in the double-scaling 
limit by using the trans-series solution to the string equation. It is then natural to ask if there 
is such a direct way of computing multi-instanton effects in the full matrix model, away from 
the critical point. Indeed, it is very easy to lift the computation in terms of the string equation 
to the original matrix integral. Recall that the 1/iV or g s expansion of the full matrix model 
is described, in the formalism of orthogonal polynomials, by a function R{z,g s ) which satisfies 
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a difference equation. In the double-scaling limit, R(z,g s ) leads to the specific heat u, and the 
difference equation satisfied by R(z,g s ) leads to the differential equation satisfied by u (the string 
equation). Difference equations, just like differential equations, also admit trans-series solutions, 
and one could suspect that the trans-series solution to the difference equation governing R(z, g s ) 
encodes the multi-instanton amplitudes of the full matrix model. To obtain the trans-series 
solutions, we consider a more general ansatz than ( |2.25| ), 

oo 

R(z,g s ) = J2C e R (e) (z,g s ), (3.19) 

£=0 



where R^°\z,g s ) is given by ( 2.25 ), and for £ > 1 we have 



RW(z,g s ) = e- £A ^/^R^(z)(l + J2^ R i,n+i(z)), i > I- (3.20) 

n=l 

Once this ansatz is plugged in the difference equation for R(z,g s ), one obtains a recursive system 
of equations for the different quantities involved. The quantity A(z), which is a parameter- 
dependent instanton action, is determined by an equation of the form 

A'(z) = /( J R ,o(^)), (3.21) 

where / is a function fixed by the difference equation. For £ = 1, n > 0, one obtains an equation 
which determines 

dRl ^ z) (3.22) 
az 

in terms of Ri n '(z) with n' < n. For n = 1, we have a differential equation for the logarithmic 
derivative, i.e. for 

The integration constant for R\^(z) can be reabsorbed in the parameter C, and for A(z) and 
the R\ >n (z), n > 1 the integration constants are fixed by using appropriate boundary conditions. 
For I > 1, the difference equation determines Ri >n m terms of Re\ n ' with £ < £' . 
In the same way, the full free energy will be given by 

oo 

F(t,g s ) = Y,C*F (l Ht,9s), (3.24) 

where 

oo 

F®{z,g s ) =e- eA ^/^ FeA (z)(l + Y,9sFi,n + i(z)), 1>1. (3.25) 

n=l 



Once ( 0.19 ) is known, one can plug it in ( 2.22 ) to deduce the F^\t,g s ). This amplitude is the 



^-instanton amplitude of the full matrix model. For example, in the case of the cubic matrix 
model, it gives log zM(t,g a ), where Z^{t ,g s ) is the partition function appearing in (pTg|) . Notice 
that, as in the case of differential equations, with the method sketched above one obtains again a 
one-parameter family of solutions parametrized by a constant C. This constant plays the same 
role as the 6 parameter £ of the original matrix model. 
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The idea of looking at trans-series solutions of the pre-string equation to obtain instanton 
corrections in the full matrix model has not been fully exploited in the literature, but it has 



been appeared in related contexts. |44|] uses essentially this approach to obtain the instanton 
action in the strongly coupled phase of the unitary matrix model. In 0], the instanton action of 
compactified c = 1 string theory is obtained by considering a trans-series ansatz for a difference 
equation of the Toda type. In the beautiful paper [78], trans-series solutions to the recursion 
equation for r n in the quartic matrix model are studied in some detail, but their focus is on the 
double-scaling limit. 

The method based on a trans-series solutions to the difference equation has two main draw- 
backs as compared to other methods. First, it does not give the value of the Stokes parameter, 
which can be computed in the saddle-point method. Second, it does not give the most general 
multi-instanton expansion for the original matrix model, since it automatically incorporates sym- 
metries of the potential. For example, in the case of the quartic matrix model with potential 
( 2.26 ) and A > there are three saddle points: a minimum at x\ = 0, and two symmetric maxima 
at X2, X3 = —X2- The most general multi-instanton amplitude will be of the form (|3.4[) 



C2 2 C3 3 Z(N 1 ,N 2 ,N 3 ), 

Ni+N 2 +N 3 =X 



(3.26) 



where we have fixed the overall normalization by setting £1 = 1. The partition function 
Z(N\, N2, N3) describes the situation where iVj eigenvalues sit at X{. However, with the method 
based on orthogonal polynomials, we find only a one-parameter family depending on a single 
constant C, 



z(c) = J2c e z 



(3.27) 



Due to the symmetry of the problem, it is easy to see that ( |3.27 ) gives the partition function 
Z((2,(3) with C2 = C3 = C, i.e. Z^ is a symmetrized instanton amplitude 



Z® = Y, Z(N-£,N 2 ,N 3 ) 

N 2 +N 3 =i 



(3.28) 



In simple cases where symmetry is not an issue, like the cubic matrix model or the GWW model, 
the method provides however the full multi-instanton amplitudes in a much more efficient way 
than alternative methods. 

It is worth noting that the relation between the expansion ( |3.19 ) and ( [3.24 ), as encoded in 



(2.22), is rather complicated. In order to extract explicit results for the free energy it is then 
useful to make it more explicit. For the £-th instanton correction, with £ > 1, it follows from 



(I2T22I ) that 



F®(z + g s ,g s ) + F®(z - g s ,g s ) - 2F®(z,g e ) 
where the superscript c denotes the connected piece, i.e. 



z,g s 



z,9s, 



S>1 



(-1) 



s-1 



£ 



R {e) (z, gs y 

RW(z,9s)_ 
R (is) (z,g s ) 



(3.29) 



i+-+e,=i 



R<f»(z,g a ) RW(z,g. 



RW(z,g s ) l^RW(z,g s )RV- k Xz,g s ) 



(3.30) 



R^(z,g s 



k=l 



(R (0 Hz,g s ))< 



+ 
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This quantity will have an expansion similar to ( p.20| 

'R {£) (z,g s 



z, 9s 



n=l 

and we obtain the following relations for the one- and two-loop contributions to (z, g t 
F^z) = \c t {z) CSC tf{^), 



(3.31) 



A'(z) 



(3.32) 



F tt2 (z) =^) + ^ooth(^) -M''(z)Qcoth 2 (^) +icsch 2 (£ 

Equations for Fg tn , n > 3 can be easily obtained from (2.22). 
3.3 An example: the quartic matrix model 

As a first example, we will study in some detail multi-instanton corrections in the quartic matrix 
model with potential ( |2.26| ) . The perturbative solution (|2.25| ) has been much studied since it was 
first worked out in the pioneering papers M, ||. The planar part is given by 



Rq,o( z ) 



A 



1 - Vl-Xz). 



As already noticed in ||, it turns out to be useful to express all results in terms of 

r = #o,q (z). 

For the higher g s corrections one finds, 

2A 2 r 



(3.33) 



(3.34) 



Roa(z) 



3 (2-Ar) 4 ' 

28A 4 r(5 + Ar) 
9 (2-Ar) 9 ' 

4A 6 r (lllA 2 r 2 + 5728Ar + 7700) 
~27~ (2 - Ar) 



(3.35) 



14 



and so on. We also recall that the double-scaling limit of R^(z,g s ) is obtained at the critical 
value A = 1 and 

9s ^0, z->l, J = {l-z)lg- 2 . (3.36) 



In this limit, 



u{K)=g s ^(2-R^\z,g s )) 



(3.37) 



satisfies the Painleve I equation ( |3.11| ) as a consequence of ( |2.27] ). 

If we now plug in the trans-series ansatz (|3.19j ) in the difference equation (2.27), we find a 
system of recursive difference equations for the R^ k \z, g s ): 

R {k) (z,9s) = ^J2 Rik ~ i] ( z >9 ~ s )( R(k) (z + g s , 9s ) + R {k) (z - g s ,g s ) + R^(z,g s )). (3.38) 
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For k = 1,2 we have, for example, 

RW(z + g s ,g s ) + RW(z-g s ,g s ) 

+ fllou^j (^Hz,g.)+RPHz + g.,g.)+RPHz-g.,g.) - y) = 0, 

R {2) (z + g s ,g s ) + R^(z-g s ,g s ) . g gg . 

+ S^K 05 ^^-) + * (0) (* + ^-) +i? (0) (^ - ft, ft) - ^ 
i?W(z,5 s )v A 

+ (* + 9s) + (* - 9s , 9s)) = 0, 

R^){z,g s )\ ) 

Using now the ansatz ( 3.2C| ) we can solve for the different quantities. Let us focus on k = 1, the 
one instanton solution. The first thing to compute is A(z), which corresponds physically to the 
instanton action. From the equation for k = 1 we find, at leading order in g s , 

e A'(z) + e - A '(z) + 4 _ I 2 . = o, (3.40) 
Ar 

which gives immediately 

caeh(A'(z)) = 2^^. (3.41) 
Ar 

This can be integrated to find A(z) up to an additive constant and an overall sign (since cosh z 
is even). Both ambiguities can be fixed by requiring that, near the critical point, 

A(z) ~ i^(l-^)!. (3.42) 



The result is 



A(z) = -j dr cosh" 1 (2^^) (l - y 

=\r(\r - 4) cosh" 1 ( — -2 J + — ^3(2 - Ar)(6 - Ar). 
4 \Ar / 2A 



(3.43) 



It can be checked that ( |3.43 ) coincides with the instanton action of the quartic matrix model 
computed in terms of its spectral curve in [36|. Notice that z stands here for the 't Hooft 
parameter. 

Once the instanton action is known, we can proceed to compute Ri t \(z). The equation one 
obtains at the next order in g s is 

R\ i (*) 1 

-±±4- = - -cotHA'(z ))A"(z), (3.44) 
Ri,i{z) 2 

which can be immediately integrated as 

R 11 (z) = [wrih(A'(z))) . (3.45) 
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The rest of the coefficients can be found by integrating the resulting equations for Ri^ n (z), and 
one finds for example, up to three loops, 

„ , . A 3 r 3 - 6A 2 r 2 + 6Ar + 24 



2^3r(2 - Ar) 5 / 2 (6 - Ar) 3 / 2 

V ; V ; (3.46) 
17A 6 r 6 - 268A 5 r 5 + 1800A 4 r 4 - 5688A 3 r 3 + 6660A 2 r 2 + 288Ar + 576 

X ^ Z ' ~ 24r 2 (Ar-6) 3 (Ar-2)5 ' 

This result can be checked by using the double-scaling limit ( |3.36| ), since with the above values 

00 5 75 

1 + £ £Ri,n+i(z) - uU(k) = 1 - ^j=k-* + (3-47) 



n=l 



which are indeed the first terms of the one-instanton trans-series solution to Painleve I 
( 3. lip (in order to compare to eq. 4.40 in [36|, one has to rescale g s — > g s /y/2). 

Using the results for Ri :U , n > 1, as well as ( [3.32 ) we find for the one-instanton, one-loop 
free energy, 

_s _ 1 ,2 i 

F$(z) = -(cosh(A'(*)) " l)" 3 (cosh(A'(z)) + l)" 3 = (3.48) 



2(3(2 -Xr)J (6-Ar)4 



while for the two and three-loop contributions we have 
5A 3 r 3 - 54A 2 r 2 + 150Ar + 24 



2^/3r(2- Ar) 5 /2(6- Ar) 3 / 2 ' 
25A 6 r 6 - 828A 5 r 5 + 10008A 4 r 4 - 50424A 3 r 3 + 89028A 2 r 2 + 7200Ar + 576 



(3.49) 



24r 2 (Ar-6) 3 (Ar-2) 5 

We can now compare these results to those obtained in f66[| . The one-loop calculation in [66] 
computes F^(z), times a coefficient. This coefficient is precisely the Stokes parameter S that 



gives the discontinuity as we cross a Stokes line. Comparing ( |3,4Sj ) with the result in |66|], we 
find complete agreement with the functional dependence on z, and we also find that in order to 
match the Stokes parameter we have to set 



gj. (3.50) 
The result at two loops in ( 3.49| ) fully agrees with the one presented in [66]. 

4. Nonperturbative effects in the unitary matrix model 
4.1 1/N expansion and phase structure 

We will first review some well-known aspects of the 1/N expansion of unitary matrix models. 



For a detailed account with many references see for example |74| . 

As in the Hermitian case, the planar limit of a unitary matrix model is described by a density 
for the eigenvalues of the unitary matrix, p(4>), which verifies the normalization condition 

p(<P)d<P = 1- (4-1) 
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The free energy of the unitary matrix model has a g s expansion of the standard form ( 2.11 ), and 
the planar free energy Fq can be computed in terms of the density as 



t 



2 r-K 



F = t d<t>p(<f>)V(4>) + - / d<t> d^p(0)p(^)lo, 



4 sin' 



(j)-ip 



(4.2) 



where t = Ng s is the 't Hooft parameter and £ is a Lagrange multiplier which imposes the 
constraint (44). The density p satisfies the equation 



^ (<^) + / dip p(V') log smZ 



- V 



+ £ = 0, 



(4.3) 



0- V 



(4.4) 



which implies that the effective potential 

V eS (4>) = -V(<P)-t ^ dV/>(V)log Mir 

J — 7T 

is constant on the support of p. 

As first found in [^7|], unitary matrix models have a rich phase structure (see [57, 32 1 for a 
detailed discussion). In the so-called ungapped phase, the density of eigenvalues has its support 
on the entire circle and is of the form (6^, 57] 



^) = ^(l + E(^ + c.c.)). 



(4.5) 



In the one-gap phase, the support of the density of eigenvalues is a single, connected interval 
inside [— w, 7r], and the density of eigenvalues is of the form 



1 a / . 2 a c . 2 a 

P{a) = 2^ff(«)cos - Wsm — - sm -, 

and is supported on the interval 

C = [-a c ,a c ] C [-it, ir]. 
The effective potential can be computed as 



V e s(4>) - V cS (a c ) = t j dzg(z) cos - y sin 2 - - sin 2 -y. 



(4.6) 



(4.7) 



(4.8) 



A simple model which exhibits two phases is the famous GWW model [[17], |8(J. The potential 
is simply 

V{z) = -{z + z~ l ). (4.9) 



In the ungapped phase, the density of eigenvalues is 



Pi<f>) 



1 

2^ 



1 

1 + j COS ( 



(4.10) 
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It is easy to see that this density is positive as long as t > 1, which is strong 't Hooft coupling. 
Therefore we will refer to the the ungapped phase as the strongly coupled phase of the model. 
For t = 1, (|4.10 ) vanishes at <fi = ir, and for t < 1 it becomes negative around (j> = tt and it is no 
longer acceptable as a solution. There must be a phase transition at t = 1, and for t < 1, i.e. in 
the weakly coupled phase, the density of eigenvalues takes the form (p~6|) with g(z) = 2/t: 



= i cos (l)V*- sill2 f' 4<1 ' (411) 

for 4> £ [-«cittc]i while it vanishes outside this interval. The endpoint of the support is deter- 
mined by the condition 

sin 2 ^ = t. (4.12) 

In Fig. |2]we show the form of the density p(<fi) as we go through the transition at t = 1. 

Using the above densities we can easily calculate the planar free energies in both phases, and 
one finds 

F w (t) = |(logt-|)+t, t<l, 

2 V l> (4.13) 

F*[t) = I, t > 1. 

Since the free energy and its two first derivatives are continuous at t = 1, we have a third order 
phase transition at large N |^. B^] . 

The GWW model can be studied as well by using the method of orthogonal polynomials 
1 44]. One first derives a recursion relation for the coefficients f n appearing in ( 2.35| ), which reads 

9s{n + l)/ n = ^(1 - / n 2 )(/n+l + fn-l)- (4.14) 

In the continuum limit 

fn^f(z,g s ) (4.15) 
and the recursion ( 4.14| ) becomes the difference equation 

(z + g s )f(z,g s ) = ^(1 - f(z,g s ))(f(z + g s ,g s ) + /(* - g s ,g s )). (4.16) 
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In order to compute the partition function W6 need the continuum limit of the coefficients r n , 



R(z,g s ), which is related to f(z,g s ) by the continuum counterpart of ( 2.36| ), 

R(z,g s ) = l- f(z-g s ,g s ). (4.17) 
One then deduces the following difference equation for R(z,g s ), 



?a/1 - R(z,g s ) = ^R(z,g 8 ){y/1 - R(z + g s ,g s ) + y/l - R{z - g s ,g s )\. 



(4.18) 



This difference equation has two different solutions depending on the value of z, which reflect 
the existence of two phases in the model [ 44 ] : 



R^\z,g s ) = \ T ^^ a{z)gf ifZ<1 (4.19) 
1 if z > 1. 



For the solution in the region z < 1 one finds 

Ro,o( z ) = z 



, (1-z) 2 ' (4.20) 
9z(z + 3) 



128(1 - z) 5 ' 



and so on. The reason for the existence of these two solutions can be easily understood if one 
looks at ( p.36[ ). This equation implies that < r n < 1. At leading order in g s , the solution 
Ro,o( z ) = z already violates this constraint, therefore when z = 1 the bound is saturated and we 
must have i?o,o( 2; ) = 1- One then finds that there are no g s corrections to this solution, and one 
ends up with (|1|). 

The method of orthogonal polynomials gives as well a very efficient way of computing the 
perturbative free energies in the weakly coupled phase. In [44] results up to genus 2 were obtained, 



but going to higher genus is just a matter of CPU time. It is convenient to compute the normalized 
free energies, where one subtracts the free energies of the Gaussian matrix model. In this way 
one obtains, 

F (t) = t, 

Fi(*) = -ilog(l-t), 

F 2 (t) - 3t ^ 



128(1 -t) 3 ' 
9t(5 + 2t) 



1024(1 -t) 6 ' 

and so on. This results are valid for t < 1. 

It is usually stated in the literature that the free energy in the strong coupling phase is given 
by its planar part, plus nonperturbative corrections coming from instantons. In fact, both phases 
have instanton corrections, albeit of a different character, and we will study both of them by 
considering trans-series solutions to the difference equation ( |4,18 ). 
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Figure 3: The effective potentiai for the GWW model, for 't Hooft parameter t = 1/4. The support of 
the eigenvalue distribution, where the effective potential is constant, is the interval [— 7r/3, 7r/3]. 



4.2 Instanton corrections in the GWW model 

We now use the method explained in section 3 to compute the instanton corrections in both 
phases of the GWW model. 

We first consider the weakly coupled phase. In this phase instanton corrections have an easy 
interpretation in terms of eigenvalue tunneling. The potential for the eigenvalues of the unitary 
matrix is —V{6) = — cos 9, which has a minimum at 9 = and a maximum at 9 = tt = —tt. The 
effective potential taking into account the eigenvalue repulsion can be computed as 



/ s i n I \ 

V cS (0) ~ V eS (a c ) = 4^-|,j (4.22) 



where 

<Z>(x) = -xy/x 2 - 1 - - cosh _1 (x). (4.23) 

This effective potential is of the form shown in Fig. |3| and its maximum is still at 9 = tt. Therefore, 
there will be multi-instanton configurations obtained by taking t eigenvalues from the support 
of the density of eigenvalues [— a c , a c ], centered around the minimum at 9 = 0, to 9 = tt = —tt, 
as shown in Fig. f|. The action of such an instanton can be easily computed, as a function of the 
't Hooft parameter, to be 

A(t) = V cS (7r) - V cS (a c ) = 4i$(H). (4.24) 
In order to find the perturbative expansion around these instantons one can use the approach 



followed in [66|. In fact, it is possible to map the unitary matrix model into a Hermitian matrix 
model [68] where the calculations of |66| could in principle be used verbatim. However, when one 
does this, the saddle point at 9 = tt is sent to infinity and this leads to extra subtleties in the 
evaluation of the saddle-point integral. On top of that, going beyond the one- instanton, two- loop 
calculation of [S6| seems rather hard. 

Instead of following this strategy, we will obtain a trans-series solution to the difference 
equation ( |4.18| ), which makes possible to calculate any instanton amplitude. By plugging the 
trans-series ansatz ( [3.19 ) in ( |4.18 ) we obtain, as before, a series of recursive relations for the 



R( k \z, g s ). For k = 1 we obtain the difference equation, 



RW(z,g,){ , y== + V 1 " R(0) ( z + 9s,9s) + A/ 1 " ^ (0) (* " 9s,9s) 
Vl - R\y)(z,g 8 ) J 

( 4 - 25 ) 

- -R^(z a ){ R(1) ( z + 9s,9s) RW(z-g s ,g s ) | = 

2 ' asi Vl-#)(, + ffs , 9s ) y/l-RW{z-gs,9sV 
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9 = 



Figure 4: An ^-instanton configuration in the weakly coupled phase of the GWW model. The eigenvalues 
of the unitary matrix live on a circle parametrized by an angle 9 G [— ir,ir], and the potential is — cos 9. 
There are N — I eigenvalues sitting around the minimum at 9 = 0, while i eigenvalues sit at the maximum 

at 9 = 7T = — 7T. 



This can be easily solved for (z, g s ), to all orders in g s . First, we solve for A(z), which is 
determined by the leading order of ( f4,25| ) as a series in g s . We find 

cosh(A'(z)) = - - 1, (4.26) 

which gives 

A(z) = -z cosh" 1 (ijz - l) + 2Vl - z, z<l. (4.27) 



As in the quartic matrix model, we have fixed the ambiguities by requiring the right behavior at 
the critical point z = 1. This agrees with the result ( |4.24| ) obtained with the effective potential 
for the eigenvalues. At order 0(g s ) we get, 

which can be immediately integrated to 

R hl (z) = (1 - sinh"5(A'(z)), (4.29) 
up to a multiplicative integration constant that can be reabsorbed in C. Since 



smh(A'{z)) = -y/(2/z -1) 2 -1 = --VT^l. (4.30) 

we finally obtain 

R 1A (z) = z^ (1 - z)^ . (4.31) 

As we explained above, we have in fact a one-parameter family of solutions parametrized by a 
constant C as in ( |3,19 ). We will determine the corresponding value of the Stokes multiplier in the 
next subsection, by comparing the instanton amplitude to known results in the double-scaling 
limit. 
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The functions Ri jVl , corresponding to n-loop amplitudes around the instanton solution, can 
be now computed in a straightforward way to any order n. We present results for n = 2, 3,4: 



3z 2 



12z - 8 



96(1 -z) 3 /V 

81z 4 - 2376z 3 + 2400z 2 + 192* + 64 



30375z 



(i 



18432(1 - zfz 2 

- 208980z 5 + 281880z 4 - 4078080z 3 + 289728z 2 - 343296z + 71168 



(4.32) 



26542080(1 - zfl 2 : 



which are valid in the weakly coupled phase with z < 1. We can then use this result to compute 
the one- instanton contribution to the free energy, by using ( 2.22j ). One finds, up to three loops, 



F (D( 



z,9s 



1^(1 
4 V 



-3/4 



1 + 



3z 2 - 60z - 8 



9s 



+ 



96(1 - zfl 2 z 

81z 4 + 79 2z 3 + 17376z 2 + 960z + 64 
18432(1 - zfz 2 



9s + 



(4.33) 



-A(z)/g a 



Z <1, 



where z stands here of course for the 't Hooft parameter and A(z) is given by (4.27). 

We now consider the strong coupling phase of the model for z > 1. This is an ungapped 
phase where the eigenvalues fill the circle, therefore we can no longer interpret instanton effects 
in terms of eigenvalue tunneling. We can however find a trans-series solution to the difference 
equations in this phase (and in fact, some ingredients of this method were already sketched in 
iH ). Since we are expanding around a different perturbative solution, the difference equations 
for the instanton amplitudes change. We obtain, for k = 1, 



4z 2 R^(z,g s ) = rW(z + g s ,g s ) + R^(z - g s ,g s ) + 2^RW(z + g s ,g s )RW(z - g s ,g s ). (4.34) 

Again, we solve first for the equation determining the instanton action, which in this case is 
simply 

(4.35) 



cosh 



and leads to 



A(z) = 2zcosh _1 (z) - 2V z 2 - 1, z > 1, 



(4.36) 



in agreement with the result of 
the result up to four loops, 



z,9s 



+ 



(z 2 - 1) 
248z 6 - 



It is also easy to compute R^'(z,g s ) to any order. We write 
2z 2 + 3 , 4z 4 + 156z 2 + 45 2 



-1/2 



1 



T9s 



12(z 2 -l) 3/2 
31716z 4 - 73602z 2 - 8505 



288 {z 2 - 1)' 



51840 (z 2 - 1) 



9/2 



9s + 



(4.37) 



Of course, by using the above result it is straightforward to compute the free energy. We also 
note that, although we have given explicit results for the one-instanton amplitude only, it is 
straightforward to calculate higher instanton corrections. 

One important aspect of our result is that the instanton action vanishes at the critical point 
z = 1 (see Fig. We can then say that the third order phase transition discovered in |47, |80H is 
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Figure 5: The instant on ac tion in the GWW model, as a function of the 't Hooft para meter t. For 
< t < 1 it is given by ( 4.27 ), and for t > 1, after the phase transition, it is given by ( [4.36 ). The action 
vanishes at the critical point 4 = 1, and it is real and positive for all t > 0. 



triggered by instantons, as first conjectured in J60|] , since instanton corrections become of order 
1 at the transition point. This instanton-driven transition has been found before in two models 
closely related to the GWW model. The first one is the quantum-mechanical unitary model 
introduced and solved by Wadia in |^] (this is the main example studied by Neuberger in fl69| ). 
The second one is continuum two-dimensional Yang-Mills theory on the sphere [fi5|] . In fact, 
it is interesting to observe that the action of the instanton in the weakly coupled phase of 2d 
Yang-Mills theory is identical to (4.27) after the identification z = A/tt 2 , where A is the area of 
the sphere and z is the 't Hooft coupling in the GWW model (we note that in order to compare 
with [ 46 1 , where the exponentiated instanton action is written like exp(— NA), one has to divide 
( 4.27| ) by z, since in our conventions we write it like exp(— A/g s )). 



4.3 Multi-instantons and the double-scaling limit 

We now discuss formal trans-series solutions in the double scaling limit of the GWW model, 
which will provide (among other things) a check of the above results. 
The double-scaling limit of the GWW model is defined by 



0, t 



1, 



9s 



2 

Ml 



t) fixed, 



(4.38) 



and describes the universal scaling near the third-order phase transition. In this limit, the 
function 

u(K)=g7^f(t,g s ), (4.39) 



where f(t,g s ) is defined by ( |4.15|) , satisfies the Painleve II equation 

u"(k) - 2u 3 {k) + 2ku{k) = (4.40) 

as a consequence of the difference equation ( |4.16| ). The double-scaled free energy F^n) is defined 
as the double-scaling limit of 

F s (t,g s )-F™(t,g s ) (4.41) 



and satisfies 



(4.42) 



Notice that the regions k, — » ±oo are mapped to t — > 1^, therefore they correspond to the weak 
and the strong coupling phase, respectively. 
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We first discuss the double-scaling limit of the weakly coupled phase. Since 

f(z,g s ) = VT^ + 0(g s ), (4.43) 

it follows that the solution to Painleve II which describes the double-scaling limit of the unitary 
matrix model must behave like 

u(k) ~ y/it, K — ► oo. (4.44) 

This asymptotic behavior determines a unique formal solution to ( [4.40 ) of the form 

(0U , n 1 73 10657 13912277 . A A . 

u w (k) = Vk g- — — H , k^oo. (4.45) 

16 ki 512 kT 8192 k~ 542888 k~ 

As in the case of Painleve I, one can consider as well exponentially suppressed corrections to this 
perturbative behavior and construct a formal trans-series solution with the structure, 



t («) = £cV>(«) = ^£c^-x e - MK3/2 eW( K ), k — ► oo, (4.46) 



1=0 e=o 



where 



and 



A = \ (4.47) 

oo 

e^( K ) = ^2u e , n+lK ~ 3n / 2 . (4.48) 

n=0 

As before, we normalize the solution with u\ t i = 1. The perturbative part u^°(k) is given by 
(4.45). The instanton expansions can be easily found by plugging the trans-series ansatz in the 
Painleve II equation. One finds a recursive equation of the form, 

{u^)" + 2ku^ -2 u^V^V* 3 ) = 0. (4.49) 

For example, satisfies the linear equation, 

(u^)" + 2/€«« - 6(u<®) 2 vP> = 0, (4.50) 
while vP> satisfies the equation 

(u^)" + 2ku^ - 6(u^) 2 u^ = 6u(°)(«W) 2 , (4.51) 
and their asymptotic expansion as k —* oo are given by 



e (D (K) = l_lZ K -3/2 + i^_ K -3 

v 1 96 18432 

e ( 2 )( K ) = i- + 5461 K -3. 
v 1 2 96 9216 



(4.52) 



( [1.46 ) gives a one-parameter family of formal solutions to the Painleve II equation which 
includes exponentially small corrections. Since the positive real axis for k is a Stokes line for 
this problem, we can ask what is the value of the Stokes parameter. This can be obtained by 
various methods. One option is a direct one-loop computation in a matrix model whose critical 
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behavior is described by Painleve II. Such a computation has been done in [|58|] for the symmetric, 
quartic matrix model with two colliding cuts. Alternatively, one can use results in the theory of 



isomonodromy deformations 55, 39]. One finds 



2tt 



(4.53) 



We can now use the results about the trans-series solution of Painleve II for k — > oo to test 
some of the results that we obained for the unitary matrix model in the weakly coupled region. 
As a consequence of ( |4.17 ), we have in this region, 



R^(z,g s ) = -2f^(z,g s )f^(z,g s ), 



(4.54) 



where f^\z,g s ) is the ^-instanton contribution to the full f(z,g s ). Therefore, in the double- 
scaling limit we should have that 



R {1) (z,g s ) -> -2it<°)(«)it (1) (K). 



(4.55) 



Indeed, one verifies from the explicit results presented above that the instanton action (|4.27|) 
behaves like, 



z -> V 



(4.56) 



and that 

(4.57) 

in agreement with ( 4.54p . Moreover, we can fix in this way the Stokes multiplier for R(z,g s ) in 
the weakly coupled phase 



Finally, one can check that 

oo 



n=l 



7T 



17 ,/2 361 

fC ' -I / 

96 18432 



(4.58) 



, 791441 q/ n 

" 3 k' 9/2 + ■■■ . (4.59) 

5308416 V ; 



Let us now discuss the double scaling limit of the strongly coupled phase. In this phase, due 
to flOg) and ( |4.17|) we have that f^(z,g s ) = and 



f^(z,g s ) = ^R^(z,g s ). 
The instanton action Q4.36 ) behaves like, 



A( z )~^(z-l)i, z->T 



therefore the relevant solution to Painleve II must behave as, 



where 



u(k) ~ e -^~ K ) 3/2 , k 



2V2 



-00, 



(4.60) 



(4.61) 



(4.62) 



(4.63) 
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As before, we can consider exponentially small corrections to the asymptotics and construct a 
one-parameter family of formal trans-series solutions with the structure, 



„(«) = = ^^^(- K )-^e-^+ 1 ) i (- K ) 3/2 eW( K ), k - -oo, (4.64) 
i=i 1=1 

where 

oo 

^( K ) = J2^,n+i(-^r 3n/2 (4.65) 

n=0 

and we normalize u\ \ = 1. It follows that u^(k) satisfies, 

(u^\k))" + 2ku^{k) =0, (4.66) 

which is, up to normalization, the Airy equation. Using the well-known asymptotics for the 
function Ai(z) (see for example pp. 101-102) we find, 

gM( K ) = l + V(--L) A / X (-^)" 3n/2 - (4-67) 



Using this we can check the result for R^\z,g s ) in the full unitary model the strongly coupled 
region, since 

(4.68) 

So far we have only considered formal solutions to the Painleve II equation, and we have obtained 
two one-parameter families of trans-series solutions characterized by their asymptotic behavior, 
namely ( |4,46| ) and ( [4.64] ) . It turns out [53| that there is a unique, actual solution to ( [4.40 ) which 



belongs to both families and has the asymptotic behaviors ( [4.44 ), ( 4.62 ). This is known as the 



Hastings-McLeod solution to Painleve II, and it defines the double-scaling limit nonperturbatively 



[20]. In the next section we will discuss the Hastings -McLeod solution and its relation to the 
formal trans-series, as well as the extension of this structure to the unitary matrix model off- 
criticality. 

4.4 Large order behavior 

An important application of instanton calculus is the determination of the large order behavior 
of perturbation theory. As in |56[], we can now use the results on the one-instanton correction 
in the weakly coupled phase of the unitary matrix model to determine the large order behavior 
of the 1/N series. More precisely, knowledge of the one-instanton contribution F^'(z,g s ) and of 
the Stokes multiplier determines a 1/g asymptotic expansion for the genus g free energies F g {z). 
The precise formula is 

F g {z) { -L T(2g + b)(-iS w )Fi,i(z) 

(4.69) 

A(z)F h2 (z) A\z)F lfi (z) 
2g + b-l (2g + b-2)(2g + b-l) 
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where b = —5/2. 

This large order formula can be tested numerically by analyzing the sequence of F g {z) for 
sufficiently large g and by removing tails with the use of Richardson transforms. The test goes 
as follows ]66|]. We first use the finite sequence 



F g (z), g = 0,l,---N (4.70) 

to construct 

Qa{z) = W^Y 9 = 0,---,N-l. (4.71) 
If flOjD holds, this sequence should have the asymptotic behavior 

*W~A^)( 1 + i i? + (?)) ' 4 ' 72 > 

as g — ► oo. We can now use Richardson transforms to extract the value of A from the k-th 
Richardson transform of the sequence {Qg} g =o,— ,N-i, which we denote by 

Qf\z), g = 0,l,--- , N-k-l. (4.73) 

The best estimate of A{z) with this sequence is then 

AW{z) = - F =^=. (4.74) 

QN-k-l\ Z ) 

The functions F\ n (z) are extracted in a similar way. For example, for F\ t \(z) and F\ 2(2) we 
consider the sequences given by 



nA(z)^ b F g (z) 

Sl ^ z) = (-igw)r( 2g + 6) ^ Fl ' l(z) ' 



2 5 / vrA(z) 2 f+ fe F 9 



2 



(4.75) 



as well as their Richardson transforms S^ k g \z). This produces the numerical estimates 

^i*) = s$-k-xW> * = 1 > 2 ( 4 - 76 ) 

for the one-loop and two-loop functions in terms of the 1/N expansion. 

In Fig. || and Fig. ^ we compare the sequences Q g , S^g, i = 1,2 and their Richardson 
transforms for k = 1,2, to the analytic results for the inverse squared instanton action, Fn(z) 
and Fx^z), respectively, for z = 1/2. The analytic results are plotted as vertical lines. Already 
from the plots we see that the agreement between the analytic prediction and the actual large 
order behavior is extremely good. Just to give some numerical examples, the estimate for the 
instanton action coming from the fifth Richardson transform, evaluated at z = 1/2, is 

A( 5 )(l/2) = 0.5328399880 (4.77) 

while the exact result is 

,4(1/2) = 0.5328399754. (4.78) 
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Figure 6: The figure in the left shows the sequence Q g defined in ( |4.71| ), as well as its Richardson 
transforms Q g k \ k = 1,2, evaluated for = 1/2 and with N = 26. The horizontal line at the top is the 



value of the inverse squared instanton action l/A 2 (l/2), where A(z) is given by (4.27). The figure on the 



right shows the sequence S\ i9 as defined in ( 4.75| ) and its Richardson transforms S^, again for z = 1/2 



and N = 26. The horizontal line is the value of i r i i i(l/2), where Fi i(z) = z 1 / 2 (l — z) 3 / 4 is the one loop 
prefactor in (4.33). 
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Figure 7: The sequence S2, g as defined in ( 4.75 ) and its Richardson transforms S^g, k = 1,2, again for 
z = 1/2 and N = 26. The horizontal line is the value of ^1^(1/2). 



Similarly, we have, for the one and two-loop estimates as compared to the exact result, 
F 1 (5 1 ) (l/2) =0.2973018513, F$(l/2) = -2.194973650, 
^11(1/2) =0.2973017788, F1M 1 / 2 ) = -2.194977300. 



(4.79) 



This analysis confirms, indeed, that the free energies F g of the GWW model in the weakly 
coupled phase diverge factorially. It is easy to verify that the instanton action ( [4.27] ) which 
controls the large order behavior is real and positive in the full phase < t < 1 (see Fig. |5|). 
Therefore, the 1/N expansion of the free energy F in the weakly coupled phase is, technically 
speaking, not Borel summable, since its Borel transform will have singularities at points of the 
form £A, i = 1, 2, ■ • • . We will discuss this issue in much more detail in the next section. 

5. Nonperturbative effects and nonperturbative definitions 

The genus expansion and its instanton corrections give only formal, asymptotic expansions of 
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the original matrix integrals. In fact, as power series in g s , these series are badly divergent. One 
can see that 

F g {z)~(2g)\, F i)n (z)~n\, i>\. (5.1) 

On the other hand, the original integrals are in many cases well-defined. In the unitary case this is 
clear, since the integration over unitary matrices is over a compact domain. In the Hermitian case, 
the general partition function ( |3.4| ) is a linear combination of convergent integrals. The question 
we will address in this section is how to recover the original, convergent matrix integrals, from the 
formal trans-series solutions. We will focus on the unitary case, since it has been comparatively 
less discussed, but we will start by analyzing the problem in the double-scaling limit, where we 
can rely on known results in the theory of exponential asymptotics and of resurgent functions. 

5.1 Formal solutions and Borel resummation 

As so often in quantum mechanics and quantum field theory, in order to give a meaning to the 
formal series we have obtained we must resum them. There are many resummation techniques 
available in the literature (see 15] for a recent review), but since the series we have to deal with 



diverge factorially, it is natural to use Borel resummation. We now briefly review some basic 
ideas of Borel resummation and of the theory of resurgent functions. 
Let 

oo 

<l>(w) = Y j a n w n (5.2) 

n=0 

be a factorially divergent series, where 

a n ~ (J3n)l (5.3) 
The Borel transform of (j), (f>(z), is defined as the series 



(5.4) 



This series defines typically a function which is analytic in a neighboorhood of the origin. If (1) 
the resulting function can be analytically continued to a neigbourhood of the positive real axis, 
and (2) the integral 

f(w)=w-VP die-*/ wl// >(t0) (5.5) 



Jo 

converges in some region of the w-plane, the series <f>(w) is said to be Borel summable in that 
region. In that case, f(w) defines a function whose asymptotics coincides with the original, 
divergent series 4>(w), and f(w) is called the Borel sum of <j>(w). 

In many cases of physical interest (like in quantum field theory and in the examples considered 
in this paper), one finds that <f>(z) can be analytically continued but it develops singularities 
(poles or branch cuts) along the real axis. This is also a typical situation in the analysis of 
irregular singular points of differential equations. Traditionally, the appearance of singularities 
on the real axis is regarded as an obstruction to Borel summability, since the integral (|5.5| ) is 
ill-defined and a prescription has to be given in order to avoid the singularities. But one can 
still use Borel resummation in order to construct well-defined quantities, and indeed this is the 
main problem addressed in the theory of resurgent functions and in the theory of exponential 
asymptotics. In the following we will rely very much on results coming from these theories. The 
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Figure 8: The paths C± avoiding the singularities of the Borel transform from above (respectively, below). 



theory of resurgent functions has been developed by Ecalle in his monumental work [31|, but for 
our purposes the results presented in |ffq, p 



75, E7| will be enough. 



In order to avoid the singularities, we will consider lateral Borel resummations. Let C± be a 
path going from to oo and avoiding the singularities of 4>{z) on the real axis from above (resp. 
below). Typically, these paths have the form shown in Fig. ||. The lateral Borel resummations 
are then defined as 

(s±ct>)(w)=w- 1 / f3 [ dte-^^litP), (5.6) 
Jc± 

provided the integral is convergent. Notice that, even if the original series has real coefficients, 
since the lateral Borel resummations are computed by integrals along paths in the complex plane, 
they lead in general to complex-valued functions. The resummations from above and from below 
are related by complex conjugation 

Hs + = s-H, (5.7) 

where 7i is the Hermitian conjugation operator (7if)(z) = f(z). Lateral resummations play 
a central role in the theory of resurgent functions, and they have been also used to resum 
nonalternating perturbative series in a variety of problems in quantum mechanics and quantum 
field theory [41, |5(|. In particular, it has been shown that in many cases the imaginary part 



which is obtained when doing these resummations has a physical interpretation. For example, 
the perturbation series for the ground state energy of the cubic oscillator is not Borel summable, 
since the Borel transform exhibits a singularity in the positive real axis. In this case, the lack 
of Borel summability is just reflecting the instability of the potential. The lateral resummations 
lead in this case to a complex answer for the Borel sum of the energy, but this is as it should be, 
since the energy of the ground state should have an imaginary part which gives the width of the 
level (see for example || [l5|] ) . 

We are interested in the divergent series which appear in the context of matrix models, 
namely when solving a differential or difference equation. Let us first focus on the case of 
differential equations. We will have in mind the cases of interest for us, namely the Painleve I 
and II differential equations. The results which we will use were obtained in [16|] and 

they hold for a large class of nonlinear differential equations with an irregular singular point at 
infinity (see [|l^] for a precise statement of the theorems and their conditions of validity) . 

As we have seen in the case of Painleve I and II, along the Stokes line arg(«) = one can 
construct a family of formal trans-series solutions of the form 



(5- 
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where 

U W(«) = K a+e Pe- eAK \W( K ), (5.9) 
a, /?, 7 are characteristic exponents of the differential equation, ^4 > is a constant, and 



(«) = ^^, n+ i/c-^ (5.10) 



e 

n=0 

are asymptotic series diverging like U£ n +i ~ n!. The Borel transforms of these series can be 
defined as above, by identifying z = k -7 , and they have singularities at the points in the positive 
real axis of the form £A, with t = 1, 2, • • • . We can avoid these singularities by performing lateral 
Borel resummations of all the formal power series appearing in the formal solution. In this way 
we construct the functions 



u ± (k;C ± ) =£c£«S?(k) = £c|^e-^( S±e W)(K). (5.11) 



It turns out that, if ReK > is sufficiently big, the infinite sum over £ is convergent in some 
angular sector around the real axis. Therefore, the lateral Borel resummations produce true, no 
longer formal, solutions of the original differential equation in this sector. Moreover, any solution 
of the differential equation with the asymptotics given by the formal series u^ ' (k) as k — > oo can 



be represented in the form (5.11) for some C± [18]. This is an important point, since we could 



think of many ways of avoiding the singularities of the Borel transform along the real axis, by 
choosing different contours. However, the general results on this type of equations tells us that 
the use of lateral resummations in the way we have explained is already enough to generate all 
the relevant solutions. In fact, we have already too many solutions, since the two contours C± in 
Fig. H produce two different families with the same asymptotic behavior, therefore they should 
be related. The relation is given by (see |L(], |i~8| ) 

u+(k;C) = U-(k;C + S), (5.12) 



where S is the Stokes parameter (which is purely imaginary). This equation gives an infinite 
number of relations between the functions u± which can be obtained by taking derivatives on 
both sides w.r.t. C and then setting C = 0. In this way one finds, 



k=l 



This expresses the fact that the lateral Borel resummations of the ^-instanton correction differ by 
an imaginary part which is exponentially suppressed with respect to their real parts. Moreover, 
at leading order this imaginary part is proportional to the (£ + l)-instanton correction. For 
£ = 0, 1 we obtain, for example, 



uf) - n (_ 0) = Su W + ■■■ = |(u« + u W ) + ■ 
n« - nL 1} = 2Su^ + ■■■ = S(uf + u^) + 



(5.14) 



where the dots denote higher order instanton contributions. 

It is clear that the one-parameter families of solutions (5.11) are in general not real, even if 



the starting point were divergent series with real coefficients. However, the resurgence relation 
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( 5.12 ) makes possible to construct a one-parameter family of solutions which are manifestly real 
for This family is given by 



u R ( K] C) =u + (k;C - S/2), C£ 



(5.15) 



To see this, notice that (5.12) gives, for C — > C — S/2, 

u+(k;C-S/2) =u-(k;C + S/2). 
therefore, if we use (|5.7|) we obtain 



(u r (k;C))* = u_(k;C + S/2) =u r (k;C), k£R, 



(5.16) 



(5.17) 



so ( |5.15| ) is real. In Ecalle's theory, the solution ( |5.15| ) is called the median resummation of 
the formal trans-series ( |5.8| ). Of course, in this solution all imaginary parts coming from Borel 
resummation cancel in the end. We have 



[k;C) = uf(K) + (C- S/2)41\k) + ■ 



so at first order in the exponential factor e Ak "' the imaginary part of (5.18|) is given by 



Im (k) + i^"Re u^) 



(5.18) 



(5.19) 



which cancels due to the first relation in (5.14). Higher order imaginary terms also cancel, and 
using these cancellations we can write the expansion of ( [5.15D , up to three instantons, in a form 
where the reality properties are manifest, 



«*(«; C) = + nL 0) ) + + ««) + ~ (C 2 - («?> + ««) 



S_ 2 

4 



(5.20) 



( ( 3 ) i (3)\ , 



The cancellation taking place in ( |5.19| ) is in fact a particularly clean example of the so- 
called cancellation of nonperturbative ambiguities. As we have seen, after Borel resummation, 
•w^) picks an imaginary part which is ambiguous and depends on the choice of contour. In the 
case of the lateral resummations we have considered, this is a sign ambiguity. This ambiguity 
is accompanied by a similar ambiguity for the coefficient of the one-instanton contribution 
which is C ± S/2. The relation ( 5.1€|) tells us how these two ambiguities should be correlated in 
such a way that they cancel in the final, real answer ( 5.15| ). 

The cancellation of nonperturbative ambiguities has been much discussed for renormalons 
1 22, |48[1 as well as for instantons in quantum mechanics (see [S3 for the most updated study and 
references to earlier work). In the quantum-mechanical double well |$3|, the standard pertur- 



bative series of the energies is not Borel summable, yet the lack of Borel summability is not a 
manifestation of an instability in the system. Therefore, there must be explicit nonperturbative 
contributions which cancel the imaginary parts incurred in when performing the Borel resum- 
mation. The ambiguous nonperturbative effect occurs in this case at the two-instanton level, 
and the final resummed answer is real. In the context of renormalon physics, the role of u(n) 
is played by a QCD observable (typically a current-current correlator where one can apply the 
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ITEP sum rules), u^°\k) corresponds to the perturbative series, and u^\k) correspond to non- 
perturbative contributions due to condensates. Consistency of the QCD path integral requires 
that the imaginary part of the Borel-resummed perturbative series cancel against the imaginary 
part of the first nontrivial condensate. In fact, this fixes the ambiguous imaginary part of the 
first condensate once a prescription is chosen for resumming the perturbative series. 

In our case, the cancellation of nonperturbative ambiguities which occurs for the median 
resummation ( |5.15 ) is a consequence of the resurgence relation ( |5.12| ) and guarantees that the 
final answer will be real. However, it does not determine the constant C. To fix C one needs 
further nonperturbative input. We will now see how the value of C can be fixed in the unitary 
matrix model and its double-scaling limit. 

5.2 The case of Painleve II 

We will now apply some of the above results to the case of Painleve II. From the trans-series 



solution (|4.46 ) we can find, using Borel resummation, two one-parameter families of solutions 
u±(k;C) to the differential equation ( 4.4CQ , and all the members of these families have the 
right asymptotics (|4.44j ) as k — > oo. The general theory sketched above tells us that these two 
families are related through the resurgence relation ( 5.12| ), which leads to the cancellation of 
nonperturbative ambiguities. 

It is instructive to verify explicitly the relations ( |5.14| ) by using resummation techniques. In 
order to do that, we compute the asymptotic expansions e^(re), up to a given order n, and we 
form the Borel transformed series 



g(0) 



(*) = £ 



9=0 



(*) = £ 



U£,n+1 



n=0 



n 



(5.21) 



z n , t>\. 



As we mentioned before, one needs to perform an analytic continuation of this series in a region 



including the positive real axis in order to be able to compute the Laplace-Borel transform (5.5) 
One way to do this in practice is to use Pade approximants. We recall that, given a series 



S{z) = ^2a k z k 



(5.22) 



fc=0 



the Pade approximant [l/m] is given by a rational function 

po +piz H \-piz l 



\l/m] S (z) 

where is fixed to 1, and one requires that 

S(z) - [l/m] s (z) = 0(z l+m+1 ). 



(5.23) 



(5.24) 



This fixes the coefficients involved in (|5.23| ). Given a series (j>(z) we can construct the Pade 
approximant of its Borel transform 



[n/2]/[(n + l)/2] 



(5.25) 
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K 


Ttt1 „(0)/ N 

lmu_|_ (^AtJ 


—ID tteuY \K>) 1 & 


1 


-0.0457932 


-0.0457633 


2 


-0.0036383676 


-0.0036383581 


3 


-0.000143729160176 


-0.000143729159748 


4 


-3.2181810964596 • 10~ 6 


-3.2181810964557 • 10" 6 


5 


-4.409270574264102 • 10" 8 


-4.409270574264109 • 10~ 8 



Table 1: Cancellation of nonperturbative ambiguities at the one-instanton level for various values of the 
double-scaled parameter k. At leading order in the trans-series expansion, lmu^\n) should be equal to 
— LSReu^ (/c)/2, as this table shows. The difference between both quantities is a three-instanton effect. 



where [•] denotes the integer part (this is the approximant proposed in [56, [15]]). V%(z) is a 
rational function with various poles on the complex plane. If the Borel transform has for example 
a branch cut, the Pade approximant will mimick this by a series of poles along the cut. The first 
pole of the approximant will be close to the branch point of the Borel transform, and increasingly 
so as n grows. A good approximation to the Borel resummed series will then be an integral of 
the form (|5.5[) where one integrates instead Vn(z). In our case, we compute 



( S± 4° ) )W = i r/ dte-^hf(t), 
U)2 Jc± 

(s±e^)(w) = - [ dte'^rfit), 
w J c+ 



VJ2 Jc± 

(5.26) 

> 1, 

c± 

where we take for C± a path from to oo along the directions ±7r/4 (like for example in |^l[ ). 
By contour deformation, these paths should give the same result as the paths shown in Fig. ||, 
at least for the true analytic continuation of the Borel transform. The Pade approximant can 
have spurious poles away from the real axis, and in some situations one might want to correct 
for these (see [p6]| ). In our case, however, all the poles of the Pade approximants are on the real 
axis, so the above integral should give a good approximation to the true result if n is sufficiently 
large. 

In table |] we compare numerically Re (k) to —iS Im uy (k)/2. All the computations have 
been done with n = 100 terms in the Borel transform, which then is used to compute the Pade 
approximant. Since we keep n fixed, our approximation will be worst as k gets small, but the 
results in the tables are exact at the level of precision that we display. We see that both results 
are quite close but do not quite agree. This is not surprising, since they are equal only up to 
three-instanton corrections which we have not calculated (these are the corrections to the first 



line in ( 5.14 )). These corrections go like exp(— 4k 3 / 2 ) and they become less important as k grows, 



as shown in the table. In table ^ we compare numerically Reu^^/t) to — iS*Imu^(K), again for 
n = 100, with similar results. 

We know that the true nonperturbative answer to the doubly-scaled unitary matrix model 
is given by the unique, real solution to Painleve II with asymptotic behaviors ( [4.44 ) and (4.62). 



This is the Hastings-McLeod solution, which we will denote by khm- This solution can be found 
by numerical integration, and we display it in Fig. 9|. On the other hand, we also know that any 
real solution to Painleve II with the asympotics ( 4.44 ) is of the form ( 5.15 ). Therefore, there 
must be a real value of C, Cum, for which 

/( K ;C hm ) = «hmW. (5.27) 
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K 


Im v$t) {k) 


—id Ke U_|_ [K ) 


1 


-0.008163 


-0.008152 


2 


-0.00004143955 


-0.00004143932 


3 


-5.53260679 • 10~ 8 


-5.53260675 • 10" 8 


4 


-2.459781001 • 10" 11 


-2.459781082 • 10" 11 


5 


-4.1875843088 • 10~ 15 


-4.1875852452 • 10" 15 



Table 2: Cancellation of nonperturbative ambiguities at the two-instanton level. 




Figure 9: The Hastings-McLeod solution to Painleve II. As k 
it decays exponentially like exp(— 2v / 2(— k) 3 / 2 /3). 



oo it asymptotes y/lz, while as . 



It turns out that this value is just Chm = 0. One way to see this is to notice that the one- 
parameter family u k (k;C) is the family of tronquee solutions studied by Boutroux In this 
family there is one single member with the asymptotics (4.62), which is the Hastings-McLeod 
solution. The value Chm = can then be obtained by comparing the structure of u r (k; C) with 
the results of [39] for the Stokes parameter of the tronquee solutions. 

We conclude that the Hastings-McLeod solution, which is the nonpertubative solution to 
the doubly-scaled unitary matrix model (and to many other models, see for example [ 
has the trans-series expansion 



0) 



1 



«hm(k) = u («; 0) = -(u + 



-S 2 i 



u<? +«<?) + 



(5.28) 



for k sufficiently large. Notice that, according to this equation, the true nonperturbative solution 
of the problem has a leading part, which is obtained by taking the real part of the Borel-resummed 
perturbative piece, and then it has instanton corrections starting with two instantons. Recall that 
these instantons are obtained by taking the double-scaling limit of the unstable configurations 
where eigenvalues of the unitary matrix sit at the maximum of the potential 6 = it. Therefore, 
these configurations do contribute to the physical answer for the partition function and have to 
be taken into account, as it was emphasized in [61]. 

It is interesting to verify the relation ( |5.28 ) numerically. In this test we really need numerical 
precision, since the second term in ( 5,28j ) comes from a two-instanton correction and is very small. 
The most delicate part is to obtain accurate results for the numerical integration of the Painleve 



II differential equation which gives whm(k). Because of this we have relied on the results of |72[. 
In table we compare the difference 



Reuf(K) 



(5.29) 
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K 


uhm(k) — Reu^(K) 


-S 2 Reuf( K )/4 


2 


0.000043768 


0.000043765 


3 


3.822644 • 10" 7 


3.822642 • 10~ 7 


4 


1.542393 • 10" 9 


1.542393 • 10~ 9 


5 


3.176 • 10" 12 


3.176 • 10" 12 



Table 3: The left column of this table shows the difference between the Hastings-McLeod solution «hm 
to Painleve II, and the real part of the Borel-resummed perturbative series Re u+ ■*(«;), for various values 
of k — k/23. This should be equal, at leading order, to the two-instanton effect — S 2 Rett+ (k)/4, which 
we show in the right column. The agreement is excellent, and the difference for small R is due to higher 
instanton corrections. 



to the two-instanton effect 



c2 

— Reuf{K), (5.30) 



and we evaluate it at different integer points n = k/2'a. The factor 2 1 / 3 has been introduced in 
order to use the results of [72|, who obtain uhm with a different normalization: 

«( K ) = -2W H h S rs (-2t K ). (5.31) 

Again, the results we display in table |3| are exact at the precision we have used (the results of 
[72] have a precision of 16 digits, and this leads to the decrease in sensitivity in the table data as 
k grows). The agreement is excellent, and the differences between the two quantities should be 
attributed to higher instanton corrections (in fact, as one can see from ( |5.20 ), these differences 



are four-instanton effects). Our numerical results confirm indeed that Chm = 0. 

To summarize, we have shown that the formal trans-series solution to Painleve II can be ap- 
propriately resummed to obtain a real, one-parameter family of true solutions with exponentially 
suppressed corrections due to multi-instantons and where nonperturbative ambiguities cancel. 
The parameter can be fixed by using further information from the nonperturbative result, and in 
particular we can reconstruct the Hastings-McLeod solution at least when k is sufficiently large. 
This gives a semi-classical expansion of this solution which includes instanton corrections. 

An interesting spinoff of our discussion is that the real part of the Borel-resummed per- 
turbative series does not give the exact nonperturbative result, since there are higher instanton 
corrections in ( |5.28| ) starting at the two-instanton level. This has been also pointed out for the 
case of the double-well potential in quantum mechanics, where the multi-instanton corrections 



play a crucial role in reconstructing the exact answer for the energies [82, 83]. In our case, this 



is intimately related to the nonlinearity of the differential equation encoding the exact answer: 
although the Borel-resummed series iJ± solve Painleve II, their sum does not. However, in sit- 
uations where the physical answer to a resummation problem is known to be real, it is often 
assumed that one should simply take the real part of the Borel-resummed perturbative series 
(see [|32], f?^]). Our example, as well as the example of the double- well, show that this is not 
necessarily the right answer, and that further nonperturbative corrections are needed. A case 
where the real part of the Borel-resummed series is known to be the exact result appears in [37], 



but the relevant quantity studied in that paper satisfies a linear differential equation, therefore 
the sum of u+ and is still a solution and higher corrections are absent. 
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5.3 Instantons and nonperturbative definition in the unitary matrix model 

The analysis of the previous section has produced a trans-series, formal solution for the free 
energy of the GWW model in the weakly coupled phase, of the form ( [3.24 ). 

oo 

T{t : g a )=glF{t,g s )=Y,C t ^ l \t,9 s ), (5.32) 



where we write 



g=o (5.33) 

F {e) (t,g s ) = g 2 s e-^fF eA (t)^(t,g s ), I > 1. 



and 



^(t,g s ) = Y / F i>n+1 (t)g^. (5.34) 

n=0 

The g s expansions involved here are factorially divergent series, therefore it is natural to consider 
as well the lateral Borel resummations 



1=0 

where 



F±(t,g s ;C ± ) = J2C±F i ± ) (t,g s ), (5.35) 



F%>(t,g a ) = g 2 s e- iA ^F £jl (t)^(t,g s ), I > 1. (5.36) 

Although our framework is slightly different from the one considered in the literature, it is known 
that some difference equations satisfy the same resurgence properties as differential equations (see 
for example [11]). It is then natural to conjecture that the resurgence relation ( |5.12 ) becomes, 
in the setting of the full matrix model, 

F + (t,g B ; C) =F_(t,g a ;C-S w ), (5.37) 

where S w is the Stokes multiplier ( |4.58| ). This leads to a cancellation of nonperturbative ambi- 
guities similar to the one considered before. At the one instanton level, we have 

^(0) _ ^(0) = S_ ( ^(i) + + . . . (5 38) 

We have tested this relation numerically by using Pade-Borel resummation of the sequences 
F g (t) (up to g = 25) and F\^ n (up to n = 15). We found very good agreement. We show some 
results for t = 1/2 and various values of ./V in table |j. Notice that in these calculations we have 
less precision since we computed fewer terms in the series, and we are not able to resolve higher 
instanton corrections. 

Assuming the resurgence relation ( 5.37 ) to hold, we can produce a one-parameter family of 
real solutions to the difference equations by taking 

F R (t,g s ;C) =F + (t,g s ;C + S w /2) = F_(t, g s ;C - S w /2), CeR. (5.39) 
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AT 

iv 


LmJ- + (t,g s ) 


-ibKeJ- + (t,g s )/2 


2 


0.000303 


0.000303 


3 


0.0000491 


0.0000491 


4 


7.440 • 10~ 6 


7.440 • 10~ 6 


5 


1.5219 • 10~ 6 


1.5219 • 10~ 6 


10 


1.4578 • 10" 9 


1.4578 • 10~ 9 



Table 4: Cancellation of nonperturbative ambiguities in the GWW model for t = 1/2 and various values 
of N. 



We now compare the resummation of the g s expansion and its multi-instanton corrections to 
the exact nonperturbative answer. Based on the analysis of the double-scaling limit, we expect 
that the resummed solution to the difference equation ( 5.39| ) with C = 



1 



^(t,g s ;0) = ^+^ 



;<w\2 



(5.40) 



is the exact, semiclassical expansion of the full nonperturbative free energy (as long as A(z)/g s is 
big enough). In order to test this relation, we have to determine the nonperturbative free energy 
of the GWW model. Since in the computation of the perturbative part we have subtracted the 
Gaussian free energy, we have 



N-l 



F np (N,g s ) = g\ log JJ hi - g 2 s F G , 



(5.41) 



i=0 



where the Gaussian free energy can be computed exactly in terms of the Barnes function, 

~G 2 (N + 1) 



N 2 

Fg = ~y log g s + lo£ 



(2k) 



(5.42) 



The product of the hi can be computed as a Toeplitz determinant |5| 



N-l 



log Y[ h i 



det 



i=0 



where 



h-i{^/g s ) 

2tt 



jfe,i=l,- JV 



In{z) 



I-n(z) 







dfl j n e-t-zcos0 

2vr 



(5.43) 



(5.44) 



are modified Bessel functions. We then find that the full nonperturbative answer for the free 
energy is 



F np (N,g s ) =g 2 J log det 4_;(l/ ffs 



W=l,- ,JV 



log 



G 2 (N + 1) 



(2k) 



log g s 



(5.45) 



which can be calculated exactly for any N,g s . According to our results, this exact function of 
N, g s has an asymptotic expansion of the form ( 5.40| ). 

As a partial verification of this statement, we have evaluated the first term in ( |5.4C| ), i.e. the 
real part of the Borel resummation of the 1/iV expansion, and compared it to the exact expression 
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0.520 



0.515 



0.510 



0.505 



Figure 10: Comparison of the exact nonperturbative result J rnp (iV, g s ) (dots) with the first term in the 
Borel resummed trans-series expansion Re J 7 ! (t, g s ) (continuous line). The horizontal axis represents N, 
the rank of the unitary matrix. The 't Hooft parameter is fixed to t = 1/2. Notice that both results are 
asymptotic, as N — + oo, to the planar limit Fo(t) = t = 1/2. 



(5.45) for various values of N, t. As in the test of ( 5.3§| ), we truncated the series at g = 25. 
The agreement is excellent and the difference between both values should come from instanton 
corrections starting at two-instantons, as we checked in detail in the double-scaled model. In 
Fig. [l(] the continuous line represents Re (t, g s ) for t = 1/2 as a function of 1 < N < 10 
(notice that, in the 1/N expansion, we can take N to be a continuous variable, equal to t/g s ). 
The dots represent the exact result ( 5.45| ) for the integers N = 1, ■ • ■ , 10. As N — » oo, both 
quantities asymptote the planar limit Fo(t) = t = 1/2. Notice that this test does not verify that 
the right value of C in ( 5.3S| ) is indeed C = 0, but this value is required in order to match the 
results in the double-scaling limit. 

The question of how accurate is the 1/N expansion in order to reproduce the full nonper- 
turbative answer has been an important one since this expansion was first formulated. In the 
particular case of the GWW model, a preliminary investigation of this issue was already per- 
formed in the seminal paper by Wadia |p0|| . We can now summarize our results, which give an 
extremely detailed answer to this question in the weakly coupled phase: 

• The 1 /N expansion of the unitary matrix model, in the weak coupling phase, is not Borel 
summable, and diverges factorially like (2g)\. Its large order behavior is governed by a 
one-instanton amplitude. 

• This 1/N expansion can however be resummed by using Borel transforms and lateral re- 
summations. Ambiguities coming from different contour prescriptions cancel against the 
ambiguity in the one-instanton contribution. 

• The exact nonperturbative answer for the model is given by the real part of the Borel- 
resummed 1/N expansion, plus an infinite series of exponentially suppresed corrections. 
These corrections can be computed as Borel-resummed instanton expansions, and they 
start at two instantons. 



The importance of instantons in the weakly coupled phase of the unitary GWW model was 
pointed out by Neuberger in |69f| , although no computational scheme was provided there to derive 
them. Their importance is closely related to the existence of a phase transition, since instanton 
contributions will become more and more relevant as we approach the transition point, where 
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their action A(t) vanishes. In particular, the trans-series expansion ( 5.40 ) will break down at 
t = 1, and in this sense we can say that the phase transition in the free energy is triggered by 



instantons which are no longer suppressed exponentially [69, f4q |. 

Although we have focused on the weakly coupled phase, it is in principle possible to obtain 
a similar convergent instanton expansion in the strongly coupled phase. 

5.4 The Hermitian case 

We have seen that, starting from the formal multi-instanton series, we can form a one-parameter 
family of solutions to the original differential or difference equations by considering the Borel 
resummations. The remaining parameter is fixed by nonperturbative input, and we did that in 
the case of Painleve II/unitary matrix model. 

In the case of Hermitian matrix models, and as we 
/ remarked in section 3, we can obtain convergent matrix 

/ C\ ^5-— integrals by suitably choosing the integration contours for 

the eigenvalues 



24 1 . Using the methods of section 
3 we can obtain a formal trans-series expression for this 
matrix integral, and by using Borel resummation we can 
recover (at least for a one-dimensional submanifold of the 
parameter space) the original, convergent matrix integral. 
Let us sketch how this would work in the case of a cubic 
matrix model with potential 




Figure 11: Three integration contours 
whe re th e cubic matrix integral defined 
by fl5.46) is convergent. 



V{z) 



1 3 

-z 
3 



(5.46) 



We can define a convergent matrix integral as in |24|] by 
choosing integration contours in the complex plane where 
the potential decreases at infinity. Three such contours are shown in Fig. [ll|. Of course, these 
contours are linearly dependent, since 

C 2 = Ci + C 3 . (5.47) 
The most general, convergent matrix integral obtained in this way is of the form ( |3.4|) 

Z= E Ci l <2 2 HNi,N 2 ). (5.48) 

N 1 +N 2 =N 

where N\ eigenvalues are integrated along the contour C±, and N2 eigenvalues are integrated 
along the contour C2 ■ As noticed in p4| , if we choose 



Ci,2 = ~ ± W 



(5.49) 



the resulting integral is real, since the contours are complex conjugate. 

On the other hand, with the procedure explained in section 3, we can obtain for the free 
energy of the cubic matrix model a trans-series expansion 



F(t,g s ) = J2C e F^(t,g s ), 



(5.50) 
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where F^\t,g s ) is the ^-instanton solution. From here, through lateral Borel resummations, we 
can construct a true, real, one parameter solution to the relevant difference equation 



F R (t, g 8 -C) = £(C - S/2) l F$ (t, g s ) 



(5.51) 



in analogy with ( 5.15| ). We have to match now the family ( |5.51 ) to the family of partition func- 
tio ns (|5.4§| ) with coefficients ( 5.4SQ . This can be done by requiring that the one-instanton term 
in ( 5.51| ) matches the one-instanton term for ( |5,48| ) as computed in [24]. This is straightforward 
and one obtains in this way 

C = i9S, (5.52) 



which is real since S is pure imaginary. We conclude that, for this choice of C, ( |5.51 ) gives a 
convergent series expansion for the logarithm of the nonperturbative answer ( 5.48| ) where the Ci,2 
are given by ( 5.49f ). As a further check of this relation, notice that, in the double scaling limit, 



(5.51) becomes the solution to Painleve I given by 



(5.53) 



If is not real, these are complex solutions, and in particular if = =pi/2 we should obtain 
the so-called triply truncated solutions of Painleve I |24j]. Using ( 5.12 ) we see that these triply 
truncated solutions correspond, respectively, to 



u + (k;0), u_(k;0), 



(5.54) 



i.e. the lateral Borel resummations above and below the real axis, which give indeed representa- 
tions of the triply truncated solutions (see for example section 5.6 of p| 



6. Nonperturbative effects in topological string theory 
6.1 General picture 

The free energy of topological string theory on a Calabi-Yau threefold X can be regarded as a 
quantum mechanical system with two different Planck constants. The worldsheet Planck constant 
is given by the square of the string length, 

K s = ll (6.1) 
For a fixed genus, the free energy F g {t) near the large radius limit has an expansion of the form 

F g (t)= ^iV s>n e-^/^. (6.2) 

n,i>0 

Here, the sum over rii, i = 1, • • • , b2(X), is a sum over topological sectors, or equivalently, over 

worldsheet instanton numbers, and n ■ t = X/i=i n '^i can ' 3e interpreted in the A model as 
the action of a worldsheet instanton with instanton numbers nj, and it depends on the Kahler 
parameters ij of X. In principle, we should expect a perturbative expansion in tut^ around 
the instanton, but the presence of worldsheet N = (2,2) supersymmetry implies that the only 
nonvanishing term in this series occurs at one-loop, and gives the Gromov-Witten invariant N g ^ n . 
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There is however a second, spacetime Planck constant which is the string coupling constant 



ftst = 9s- (6.3) 
Indeed, the (perturbative) free energy is given by a series of the form 

oo 

F^(t,g s ) = J2F g (t)9 2 s 9 - 2 - (6.4) 

9=0 

It is then natural to conjecture that the full free energy of topological string theory should be in 
general a trans-series expansion depending on two small parameters, namely 

9s, e- A ®/9° (6.5) 

where A(t) is the action of a spacetime instanton, i.e. 

oo 

F(t,g,) = Y i FV>(t,g a ), F® (t, g s ) ~ e ~ eA ^ , t > 1. (6.6) 
e=o 

Here we have assumed for simplicity that all spacetime instantons are classified by a single 
instanton number £, but of course there can be more general situations. 

This conjecture was put forward in |36| by using the connection between topological strings 
and matrix models, and some simple examples of topological string theory were studied there 
from this point of view. These examples are topological string theory on the toric Calabi-Yau 

X p = O(-p) 0(p - 2) -> P 1 (6.7) 

and its p —* oo limit, which can be interpreted as a theory of simple Hurwitz numbers. It was 



proposed in [S5| that these topological string models models can be described in terms of spectral 
curves akin to those appearing in matrix models, and this was shown to be the case in ||34|| . This 
makes possible to calculate explicitly the spacetime instanton action A(t), as well as the first 
few terms of the one-instanton contribution to the free energy F^ (t, g s ), by using saddle-point, 
matrix model techniques. We will now provide further evidence for this general conjecture by 
analyzing simple topological string models with the techniques and ideas developed above. 

6.2 A toy model 

In the previous sections we have shown that multi-instanton series can be obtained, in the case 
of matrix models, by finding trans-series solutions to the difference equations that describe the 
model. In general, it is not known if the free energies of topological string models are described 



by differential or difference equations. Here we point out that the Hurwitz model studied in [66] 
can be described by a difference equation [^] which admits a trans-series solution. Therefore, 
the total free energy of the Hurwitz model is indeed of the form (|6.6| ). We will also verify that 
this solution reproduces the one-instanton effects computed in J66[. 



The Hurwitz model is defined, in down-to-earth terms, by a partition function of the form 

(2g - 2 + 2d) 



g>0 d>0 
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where Q = e tH and gn can be regarded as formal parameters keeping track of the degree and the 
genus, respectively, and d (l d ) is a simple Hurwitz number counting degree d covering maps 
of P , with simple branch points only , and by Riemann surfaces of genus g (see [66] for explicit 
expressions). The free energy log Z describes connected, simple Hurwitz numbers d (l d )', 



(2a- 2 + 2d)! 



and it has the genus expansion 



F(g H ,t H ) = Y,9h~ 2f b(Q)- ( 6 - 10 ) 

3=0 

This theory is in fact a topological string theory in disguise. It can be realized as a special 
limit of topological string theory on certain toric Calabi-Yau manifolds, see for example [fl4"| , 17] 
for detailed derivations. It was conjectured in |65|] and proved in |34| ] that Hurwitz theory can 
be described in terms of matrix integrals, and this in turn was used in |36| to compute the one- 
instanton contribution to the perturbative free energy ( |6.10| ). In particular, it was found in [66] 
that the instanton action is given by 



A(t H ) =2w(x + cosh(w)x* -2), (6.11) 
In this equation, the dependence on tn occurs through the variable x defined by 

X e" x = e-* H , (6.12) 

and w is defined by the implicit equation 

w 1 . 

X*. (6.13) 



sinh(u; 



As shown in TC], the free energy of Hurwitz theory satisfies a difference equation of the Toda 

type, 

exp(F(t H + g H ) + F(t H - g H ) - 2F{t H )) = g 2 H ^ d^F(t H ,g H ). (6.14) 

As we did in section 3, we can try to solve this equation with a trans-series ansatz of the 
form ( 3.24 ). Doing this one immediately obtains the following equation for the one-instanton 
amplitude, 

e W [A gH F^(t H ))A gH F^(t H ) = g 2 H e tH dlFU, (6.15) 

where we have written 

A h /(t) = f(t + h)+ f(t -h)- 2f(t) (6.16) 

to denote the discrete Laplace operator with step h. The first term in the expansion of fl6l5|) m 
powers of gu gives an equation for A'(tH), 

co S h(A'(t H )) - l] = ^-^{A\t H ))\ (6.17) 

where we used that 

d 2 tH F^\t H ) = X . (6.18) 
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One can check that the function A(tj{) defined implicitly by ( |6.17 ) coincides with the instanton 
action (6.11) computed by the matrix model (we tested this by expanding both quantities around 
the critical point of the model at % = 1). Of course, it is straightforward to use ( |6.14 ) to derive 
a full trans-series solution for the free energy. By the arguments already explained above, one 
can perform if needed Borel resummations to obtain a one-parameter family of true solutions. 
We expect from the general arguments explained in section 4 that all solutions to the Toda- 
like equation ( |6. 14 ) with the asymptotics fixed by the perturbative expansion are described by 
this one-parameter family. Therefore, a nonperturbative completion of the theory should be 
equivalent to fixing a value for this parameter. 

6.3 Holographic description and nonperturbative effects 

In the example considered above, as well as in the more general example of topological string 
theory on local curves studied in [66], one can compute nonperturbative effects by using a matrix 
model dual description. On the other hand, there is strong evidence |65|, |l0| that the closed and 
open amplitudes of topological string theory on a toric Calabi-Yau threefold can be described in 
terms of recursion relations on a spectral curve typical from matrix models [BE|. It is then natural 
to expect that, in the same way that the full matrix model partition function involves a trans- 
series expansion obtained by summing over all instanton sectors (i.e. over all filling fractions), as 



in (3.4), the full partition function for topological string theory on a toric Calabi-Yau threefold 
will involve such a sum over spacetime instanton sectors, as we have conjectured above. 

A partial verification of this expectation, beyond the simple toy model considered before, 
comes from looking at topological strings with large N Chern-Simons duals. In particular, 
topological string theory on A p _\ fibrations over P 1 is conjectured to be dual to Chern-Simons 
theory on the lens space S 3 /Z p H, and some detailed evidence for this was obtained in 51, 50]. 
In this case, near the orbifold point in moduli space which is dual to the Gaussian point of Chern- 
Simons theory, the matrix model realization can be made explicit by using the matrix integral 
representation of Chern-Simons theory found in Q . One finds that the perturbative topological 
string partition function, near the orbifold point, is given by a matrix integral similar to (|3,3|) 



Z(ti,... ,t p ) = Z(N 1 
1 



,N P ) 



N 



In this equation, 



t— r dAj -i— r / , Xi-XjV __2^" (A— A*) 2 ( 6 - 19 ) 

n~jr ni 2sinh ~ * e 



i<j 




(6.20) 

the contour passes through the point 2ni(k — l)/p, and the Kahler parameters of the Calabi- 
Yau ti are identified with the 't Hooft parameters g s Ni. The nonperturbative answer, which is 
the Chern-Simons partition function on §> 3 /Z p , is given by [63. [jj 



Z C s(N,g. 



(^■■■Cp p Z(N 1 ,... ,N P ), 



(6.21) 



where 



Q =exp(— (i-1) 2 ), j = l, 



(6.22) 
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and k = k + N is the shifted coupling constant of Chern-Simons theory, which is related to the 
string coupling constant by 

9s = -p. (6.23) 
pk 

Notice that, in Chern-Simons theory, k is an integer. 

The expression ( |6.21| ) is precisely of the form (|3,4|). Therefore, in topological string the- 
ories on toric Calabi-Yau threefolds with large N duals, the nonperturbative definition of the 
topological string partition function, which can be read from the Chern-Simons gauge theory 
dual, involves a sum over all the instanton sectors of the matrix model ( 6.1S| ). This provides a 



further confirmation of the conjectural structure of the full topological string partition function 
as involving a sum over spacetime instanton sectors. 



In [33], Eynard has pointed out that the sum over multi-instantons in a matrix model is 
independent on the choice of filling fractions, and therefore should be background independent. 
Indeed, since in an expression like and (^2^) we sum over all possible backgrounds ti = g s Ni, 
the final result should not depend on any particular choice of background. This is in contrast to 
the perturbative topological string free energy, where one has chosen a fixed, arbitrary background 
given by ti = g s Ni. Notice that, albeit the ti transform in a nontrivial way under the symplectic 
group acting on the periods of the Calabi-Yau, the total 't Hooft coupling t = t\ + ■ ■ ■ + t p , which 



is the variable appearing in the l.h.s. of (6.21), should be modular invariant for this picture 



to be consistent. For p = 2 this can be checked by using the results of Jy]. As we have just 
seen, holographic duals force us to consider precisely the sum over instanton sectors as a natural 
nonperturbative definition of the full topological string partition function. This suggests that 
( 6.21[ ) is a natural starting point to construct background independent topological string models. 



7. Conclusions and open problems 

In this paper we have studied various aspects of nonperturbative effects in matrix models. First of 
all, we have developed techniques to compute multi-instanton amplitudes by finding trans-series 
solutions to the relevant difference equations, and studied in some detail both the Hermitian, 
quartic matrix model, and the unitary GWW model. These techniques give formal, asymptotic 
series, and by using results from the theory of resurgent functions and of exponential asymptotics, 
we spelled out in detail how to obtain convergent series which can then be used to find multi- 
instanton expansions of the true, nonperturbative matrix integrals. We illustrated this in the 
case of the unitary matrix model and the GWW model, and we clarified in this way some subtle 
aspects of the 1/N expansion which might be relevant in more complicated situations. Finally, 
following ]6q| , we argued that these trans-series instanton expansions should be also relevant in 
topological string theory, and we gave some pieces of evidence for this. In particular, we showed 
that in topological strings with both a matrix description and a holographic Chern-Simons dual, 
the nonperturbative definition in terms of the gauge theory partition function indeed forces us 
to consider all the instanton sectors of the matrix model. 

There are many aspects of the paper that should be further clarified and extended. We end 
with a list of open problems which we find interesting. 

• The strategy followed here to study Painleve II can be also used to study nonperturbative 
effects in the (p, q) minimal string. In particular, the non-unitary models (like the (2, 5) 
model that describes the Yang-Lee singularity) are well-defined nonperturbatively |l3j and 
one could "unfold" the semi-classical content of the exact answer by using the approach 
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based on trans-series solutions. Another closely related example is the weak coupling phase 
of the minimal superstring with flux, which is described by a close cousin of Painleve II 



5E , 7€]. The only nontrivial information, namely the value of the Stokes parameter as a 



function of the flux, can be inferred from the results in |3 



The instanton effects computed for the weakly coupled phase of the GWW model should be 
inherited in the multi-trace unitary models used to describe (super) Yang-Mills theory on 
S 3 x S 1 . These effects are of order 0(e N ) and they might provide tractable gauge theory 
duals to D-brane effects on the string side. 

Another model that one could study with these techniques is the nonperturbative comple- 
tion of 2d gravity proposed in |2l]| , which is described by a different string equation yet 
has the same asymptotics than Painleve I. Since the perturbative asymptotics determines 
the Stokes parameter and the one-instanton amplitude, the formal trans-series extension 



of the solution proposed in |21| must share many properties with the trans-series solution 



to Painleve I, and it would be interesting to understand their relation in more detail. 

Although the results we use for Painleve II are corollaries of more general results in the 
theory or resurgent functions and of exponential asymptotics, the results we have presented 
for the full matrix model have not been established rigorously. It would be interesting to 
show, by extending known results and techniques, that the properties we have assumed and 
tested numerically indeed hold for the difference equations characterizing matrix models. 
This is potentially a very rich arena, since all the relevant quantities in the trans-series 
asymptotics, like the instanton action, depend now on parameters (the 't Hooft coupling 
and the coupling constants of the model), and we will have a very rich situation in which 
the analyticity structure (for example, the location of the poles of the Borel transform) 
changes as we move in parameter space. 

It seems very likely that the topological string theory on the Calabi-Yau threefold X p 
defined in ( |6.7D is also described by a difference equation which generalizes ( 6.14| ). It would 



be interesting to find such an equation and use it to obtain trans-series solutions. Of course, 
it would be even more interesting to find explicit difference equations for the topological 
string partition function on other Calabi-Yau targets, or to translate the matrix model 
results of |35|, IC] in such a framework. 



As we mentioned in section 3, general instanton amplitudes are closely related to multicut 
amplitudes, and further clarification of this relationship should be beneficial in the study 
of nonperturbative effects in matrix models. Results in this direction will appear in 



As we explained in the last section, for topological strings with matrix model as well 
as holographic duals, the idea of completing the perturbative topological string partition 
function by adding instanton sectors of the matrix model is not only reasonable; it is in 
fact imposed to us by the holographic dual. It would be very important to clarify these 
nonperturbative effects and to understand their implications for background independence 



in string theory, as pointed out in 33], and more generally for large N dualities as a whole. 
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